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1. 


INTRODUCTION 


1.1  MOTIVATION  FOR  MODERN  ESTIMATION  AND  CONTROL 
IN  ADAPTIVE  OPTICS 

Adaptive  optics  concepts  and  technology  have  stimu- 
lated efforts  among  members  of  the  laser,  surveillance,  and 
astronomy  communities.  Classical  optical  systems  designed 
for  such  applications  have  reflected  the  premium  put  on  opti- 
cal precision.  As  a consequence,  these  systems  have  character- 
istically been  heavy,  larger  than  necessary,  and  very  expensive. 
They  also  have  been  incapable  of  compensating  for  the  effects 
of  the  atmosphere.  Adaptive  optics  have  the  potential  to 
ameliorate  most  of  these  drawbacks  because  they  are  designed 
to  sense  and  correct  optical  system  phasefront  aberrations. 

Such  devices  permit  relaxation  of  requirements  for  optical 
precision,  thereby  reducing  their  weight,  cost,  and  size  and 
can  eliminate  blur  and  power  dissipation  caused  by  the  atmos- 
phere. 


The  method  for  phasefront  sensing  varies  among  adap- 
tive optics  systems  which  have  been  implemented.  Direct  phase- 
front  measurement  of  either  backscattered  energy  from  a target 
or  target  thermal  radiation  is  possible,  as  is  the  inference 
of  phasefront  via  observation  of  the  effect  of  probe  signals 
on  signal  strength  at  a target  (active  illumination)  or  the 
strength  of  a signal  recovered  from  a passively  radiating 
source.  For  both  phasefront  sensing  schemes,  however,  the 
method  of  correcting  phasefront  is  based  upon  radiation 
received  from  the  target;  dynamic  optics  (tilt  mirrors  for 
tilt  errors,  relative  displacemert  of  secondary  and  primary 


mirrors  in  beam  expanders  for  focus  errors,  and  deformable 
mirrors  for  either  or  both  tilt  and  focus  plus  higher-order 
aberration  correction)  are  positioned  to  implement  the  cor- 
rections. For  the  direct  measurement  (phase  conjugate) 
system,  however,  the  optical  path  over  which  the  phase  mea- 
surement is  .nade  must  be  the  same  as  that  over  which  the 
correction  is  implemented.  It  is  then  assumed  that  the  con- 
jugate of  the  measurement  will  correct  for  the  aberrations. 

In  the  direct  or  intensity  measurement/probe  signal  (active 
multidither)  scheme,  the  intensity  can  be  collected  over  a 
different  optical  path. 

Possible  adaptive  optics  implementations  in  a con- 
ceptual system  for  maximizing  high  energy  laser  (HEL)  flux 
density  are  shown  in  Figs.  1.1-1  and  1.1-2.  The  phase  con- 
jugate implementation  (Fig.  1.1-1)  is  for  a passively  radiating 
target;  energy  is  collected  at  a different  wavelength  than 
that  transmitted  by  the  HEL.  The  target  radiation  propagates 
through  the  shared  optical  train  to  a beamsplitter  which  re- 
flects the  target  energy  in  a different  direction  than  the 
outgoing  HEL  wavefront.  The  split  target  radiation  is  then 
imaged,  spatially  sampled,  and  processed  to  provide  phase- 
front  corrections  which  are  implemented  by  the  dynamic  optical 
components,  and  presumably  compensate  for  all  aberrations  intro- 
duced between  the  target  and  the  sensor. 

The  phase  conjugate  concept  will  not  correct  for 
flux  density  loss  at  the  target  due  to  HEL  beam  aberrations. 

An  additional  beam  cleanup  system  is  included  in  Fig.  1.1-1 
for  this  purpose.  It  includes  a sensor  to  measure  the  HEL 
output  phasefront  and,  possibly,  a separate  deformable  mirror 
to  correct  it,  although  it  appears  that  such  corrections  could 
be  combined  with  those  from  the  phase  conjugate  device  and  im- 
plemented by  its  dynamic  optics. 


o 
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In  Fig.  1.3-2  an  active  multidither  scheme  has  been 
substituted  for  both  the  phase  conjugate  and  beam  cleanup 
systems.  Active  multidither  first  modulates  the  outgoing 
HEL  phasefront  with  probe  (dither)  signals  via  the  active 
optics.  It  then  collects  HEL  energy  reflected  from  the  tar- 
get and  determines  the  effect  of  the  dither  on  that  intensity. 
If,  for  example,  increasing  the  amplitude  (during  the  positive 
part  of  the  dither  cycle)  of  a mode  shape  implemented  on  the 
deformable  mirror  increases  the  strengtn  of  the  target  return 
slightly,  then  the  amplitude  of  that  mode  should  be  made  greater 
to  compensate  for  phasefront  aberrations.  In  this  way,  the 
active  multidither  "hill  climbs"  its  way  to  maximizing  inten- 
sity on  target  or,  equivalently,  to  eliminating  the  effect 
of  phasefront  aberrations  existing  in  the  HEL  beam  and  created 
in  the  propagation  path  between  the  beam  and  the  target. 

Adaptive  optics  concepts,  although  promising,  are  far 
from  being  devoid  of  problems.  In  fact,  the  noisy,  coupled- 
multivariable,  time  varying,  potentially  nonlinear,  high  or- 
der, and  distributed  (space/time)  characteristics  of  such 
devices  challenge  the  most  competent  controls  engineer. 


Phase  conjugate  sensors  are  noisy  and  their  noise  properties 
correlate  with  the  intensity  of  the  received  radiation,  which 
changes  because  of  KEL  beating  and  the  amount  of  aberration 
introduced  by  the  system  optics  and  atmosphere.  Ambient  radi- 
ation also  introduces  noise  in  these  sensors.  Beam  cleanup 
sensors  must  contend  with  their  internal  noise  and  the  rapidly 


varying  intensity  from  the  HEL.  For  both  sensors  the  number 
of  phasefront  samples  across  the  beam  may  be  limited  since, 
for  a given  strength  of  the  sampled  beam,  the  phasefront  de- 
tector outputs  decrease  as  the  beam  size  is  increased  and/or 
the  detector  size  is  decreased.  Active  multidither  sensor 
outputs  are  corrupted  from  internal  noise,  collected  ambient 


- i radiation,  HEL  intensity  variations  in  the  target  backscattered 
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energy  and  speckle.  Speckle  is  an  intensity  variation  char- 
acteristic of  coherent  radiation  introduced  by  the  non- 
specular  nature  of  the  target,  and  it  can  introduce  severe 
levels  of  return  signal  modulation.  Such  modulation  is  a 
function  of  target  parameters  and  beam  spot  size  at  the  tar- 
get . 


In  addition  to  the  estimation  problem  resulting  from 
such  noise  effects,  the  control  of  phasefront  aberrations  is 
difficult  because  multiple  controlled  optical  elements  (and 
potentially  multiple  coupled  controlled  elements  within  them, 
e.g.,  the  deformable  mirror)  are  required,  which  may  have 
nonlinear  dynamics.  The  actual  position  of  a deformable  mirror 
surface,  in  response  to  input  commands,  may  not  be  known.  Also, 
the  control  system  gain  in  an  active  multidither  sys=tem  can 
vary  depending  on  the  degree  of  convergence  of  the  HEL  optical 
system.  Obviously,  any  estimation  and  control  structures  for 
adaptive  optics  must  provide  good  performance  in  the  presence 
of  the  varying  noise  and  system  dynamic  properties. 


Modern  estimation  and  control  techniques  are  well 
suited  to  the  solution  of  the  problem  areas  described  for 
adaptive  optics  systems.  Kalman-type  estimators,  e.g.,  Kalman, 
extended  Kalman,  and  adaptive  filters,  can  provide  optimal  or 
near-optimal  estimates  of  phasefronts.  Correspondingly,  opti- 
mal control  methodology  can  provide  control  structures  which 
derive  the  best  possible  performance  from  a system  and 
equally  as  important,  can  define  which  variables  and  parame- 
ters are  most  important  to  providing  good  system  response 
and  how  they  must  be  varied  to  reduce  the  phasefront  errors. 


In  the  present  study,  the  three  adaptive  optics  sys- 
tem types  (phase  conjugate,  beam  cleanup,  and  active  multi- 
dither) have  been  modeled  and  optimal  estimation  and  control 
structures  have  been  designed  and  evaluated  for  each. 
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ESTIMATION  AND  CONTROL  APPROACH  ' 
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Figure  1.2-1  illustrates  the  modern  estimation  and 
control  approach  of  the  present  study.  The  environmental 
disturbances  are  characterized  by  a truth  model.  For  example, 
turbulence  spectra  provided  by  AFWL  (Ref.  1)  were  approximated 
by  linear  Gauss-Markov  models,  i.e.,  linear  filters  driven  by 
white  noise.  Those  truth  model  states  which  degrade  the  per- 
formance (target  intensity)  of  the  system  are  extracted  by 
T 

the  matrix  and  added  to  the  active  optics  outputs  which 
attempt  to  counteract  these  disturbances.  A sensor  measures 
the  disturbed  wavefront,  or  the  effect  of  the  wavefront  dis- 
turbances on  target  intensity.  A Kalman-type  filter  uses  the 
measurement,  £•,  and  the  assumed  truth  model  dynamics  to  derive 
a statistically  optimal  estimate  of  the  truth  model  states,  z. 
Statistically  optimal  implies  that  the  Kalman  filter  establishes 
the  best  balance  between  the  uncertainty  of  the  measurement  due 
to  sensor  noise  and  the  uncertainty  of  the  random  process  it 
is  attempting  to  estimate.  It  is  also  significant  that  the 
Kalman  filter  computes  an  estimate  of  the  full  disturbance 
state;  for  example,  if  a turbulence  aberration  is  characterized 
by  a third  order  Markov  model,  the  Kalman  filter  estimates  the 
magnitude  of  the  aberration  as  well  as  two  states  related  to 
the  time  derivatives  of  the  aberration.  The  derivative  states 
can  be  used  to  advantage  in  computing  controls  since  they  give 

some  indication  of  the  future  evolution  of  the  disturbance 

T 

process.  The  control  gain  matrix,  K , computes  the  optimal 
linear  combination  of  disturbance  magnitude  and  derivative 
states  to  produce  the  control  commands  to  the  active  optics 
components . 
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In  a system  design  that  involves  linear  plant  and 


disturbance  dynamics,  a quadratic  performance  index,  and 
Gaussian  noise  statistics,  the  separation  theorem  (Ref.  2) 


Figure  1.2-1 


Optimal  Estimation  - Control  Block 
Diagram  for  Adaptive  Optics  System 


allows  the  estimation  and  control  syntheses  to  be  performed 
independently.  Figure  1.2-2  illustrates  the  application  of 
this  theorem  to  the  present  'study.  In  estimation  design  and 
development,  perfect  control  is  assumed  (i.e.,  instantaneous 
response  of  the  active  optics);  the  control  design  is  carried 
out  asstiming  perfect  state  estimation  (i.e.,  estimated  states 
are  identical  to  truth  states).  In  the  context  of  modern 
estimation  and  control  synthesis,  AMD  and  PC  are  alternative 
estimation  methods;  the  optimal  control  design  approach  ap- 
plies to  both  the  AMD  and  PC  systems. 

It  should  be  noted  that  the  nonlinear  measurement 
involved  in  the  AMD  system  precludes  dogmatic  application 
of  the  separation  theorem.  Separation  of  the  estimation 
and  control  syntheses  is  a practical,  though  approximate, 
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ESTIMATION  DESIGN  SYSTEM 


CONTSOI  DESIGN  SYSTEM 


• «)  • ESTIMATION  ERROR 


5C  ■ <i  ’ CONTROL  ERROR 


Figure  1.2-2  Separation  of  Estimation  and 

Control  Syntheses 


design  technique;  computer  simulation  of  the  integrated  esti- 
mation and  control  system  validates  the  success  (or  failure) 
of  the  system. 


1.3  OVERVIEW  OF  THE  REPORT 


Chapter  2 of  this  report  outlines  the  active  optics 
and  environmental  disturbance  models  used  in  the  present 
study;  relevant  parameter  values  are  presented  in  Appendix  A. 
The  AMD  and  PC  estimation  algorithms  are  described  in  Chapter 
3;  alternative  classical  estimation  techniques  for  PC  are  also 
described.  Optimal  disturbance  nulling  control  and  its  appli- 
es cation  to  the  tilt,  focus,  and  deformable  mirror  control  chan- 

| nel3  are  covered  in  Chapter  4.  Chapters  3 and  4 also  include 

j * preliminary  performance  evaluations  for  the  estimation  with 

t perfect  control  and  control  with  perfect  estimation  cases, 

■ ' respectively.  Computer  simulation  results  for  the  integrated 

' > estimation  and  control  system  evaluation  are  presented  in 
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Chapter  5.  The  beam  cleanup  (BC)  system  is  addressed  in 
Chapter  6;  optimization  of  a three  spectral  line  control 
system  is  described,  and  simulation  results  for  a candidate 
system  are  presented.  Chapter  7 presents  the  conclusions  and 
a discussion  of  areas  for  further  refinement  of  the  adaptive 
optics  system  and  controller. 


2. 


MODELS  FOR  ADAPTIVE  OPTICS  CONCEPTS 


In  this  chapter,  the  mathematical  models  used  in 
adaptive  optics  estimation  and  control  algorithm  development 
are  described.  The  models  are  divided  into  two  major  subsets 
according  to  whether  a particular  model  simulates  an  adaptive 
optics  system  component  (such  as  a mirror)  or  an  environmental 
effect  (such  as  beam  diffraction  effects  in  propagation).  The 
majority  of  the  models  are  common  to  both  the  AMD  and  PC  systems; 
the  only  significant  differences  between  the  two  systems  are 
the  type  of  sensor  used,  the  distribution  of  thermal  distor- 
tion effects,  and  the  target  effects  modeled. 

The  models  used  in  this  study  were  simplified  to  be 
compatible  with  control  synthesis;  the  models  do,  however,  in- 
clude all  important  phenomena  observed  in  breadboard  tests  or 
implied  by  recent  analysis  results.  Parameter  values  and  con- 
trol specifications  (which  are  listed  in  Appendix  A)  have  been 
chosen  to  reflect  the  current  state-of-the-art  in  hardware 
development . 


2.1  SYSTEM  MODELS 

2.1.1  Active  Multidither 

The  component  models  and  component  interactions  for 
Active  Multidither  (AMD)  simulation  are  shown  in  Fig.  2.1-1. 
The  adaptive  optics  system  consists  of  the  high  energy  laser 
(HEL),  a deformable  mirror,  a two  axis  tilt  mirror,  a focus 
mirror,  and  a photodetector  sensing  element.  Disturbances 
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Figure  2.1-1  Active  Multidither  Model 


within  the  adaptive  optics  system  due  to  mechanical  vibration 
and  mirror  heating  are  generated  by  the  Vibration  and  Thermal 
Distortion  models,  respectively.  Atmospheric  turbulence  and 
beam  propagation  effects  are  considered  in  the  Turbulence 
and  Target  Illumination  models.  Effects  of  target  surface 


> 

J 

f 


! 


properties  are  provided  by  the  Speckle  model;  the  Photodetector 
model  includes  the  effects  of  background  radiation  and  internal 
noise.  Transit  delays  from  beam  propagation  to  the  target  and 
return  to  the  receiver  are  also  included  in  the  AMD  system 
model . 


The  phase  aberration  function,  which  represents  the 
perturbation  of  the  outgoing  wavefront  from  that  of  an  ideal 
converging  spherical  wave,  is  expressed  in  terms  of  Zernike 
circle  functions  (Ref.  1),  i.e. 


where : 


<P(x,y ) 
P^x.y) 

ci 


$U,y)  - & l c .Pt  (x,y)  (2.1-1) 

A i-1  1 1 

the  phase  aberration  function 

the  is  Zernike  circle  function 
t h 

the  i Zernike  coefficient 


The  expansion  of  Eq . 2.1-1  is  truncated  at  five  terms;  the 
P^ ' s for  these  five  terms  are  listed  in  Table  2.1-1, 


TABLE  2.1-1 

IMPLEMENTED  OPTICAL  MODES 
___ T-1755 


RELATIONSHIP 

OPTICAL  MODE  SHAPE 

pi  ■ 31 

P2  - 7 

Py  * x2  + y2 

P4  - x2  - y2  ) 

> 

P5  - xy  ) 

Pitch  Tilt 

Yaw  Tilt 

Focus 

Astigmatism 
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The  variables  x and  y are  distances  from  the  beam  center 
nondimensionalized  by  the  aperture  radius,  i.e. 

x - (2.1-2) 

aperture 

The  intensity  measurement  of  the  receiver,  the 
position  transducer  outputs  of  the  tilt  and  focus  mirrors , 
and  the  input  command  to  the  deformable  mirror  are  supplied 
to  a digital  extended  Kalman  filter.  The  state  estimates 
generated  by  the  Kalman  filter  are  the  inputs  to  the  active 
optics  components;  the  sample  and  hold  and  analog  controller 
functions  are  included  within  the  deformable  mirror,  tilt 
mirror,  and  focus  mirror  modules. 

2.1.2  Phase  Conjugate 

The  models  and  interactions  of  the  Phase. Con jugate 
(PC)  system  are  shown  in  Fig.  2.1-2.  The  component  and  en- 
vironment models  used  in  the  Phase  Conjugate  system  model 
are  quite  similar  to  those  of  the  AMD  system;  the  only  signi- 
ficant differences  between  the  two  system  models  are 

• The  receiver  model  is  a "black  box"  represen- 
tation of  a shearing  interferometer  as  opposed 
to  a photodetector  in  the  AMD  system 

• Thermal  distortion  effects  are  simulated  by 
two  distinct  models;  a separate  model  is  pro- 
vided for  those  optical  components  which  are 
part  of  the  shared  aperture  optical  train 

• Target  effects  (speckle)  do  not  influence 
receiver  performance. 

In  phase  conjugate  adapative  optics  a shearing 
interferometer  measures  the  incoming  wave  phasefront  at 
? ;j  the  shared  aperture  exit  pupil.  As  indicated  by  Fig.  2.1-2 
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Figure  2.1-2  Phase  Conjugate  System 

this  incoming  wave  is  generated  by  a point  source  model  of 
the  target;  in  addition  to  the  passive  target  radiation, 

1 aster  illumination  effects  can  influence  (increase)  the 
signal  level  at  the  receiver.  Since  the  incoming  wave 
passes  through  the  shared  apeture  optical  train  before 


reaching  the  shearing  interferometer,  it  is  affected  by 
thermal  distortions  of  the  shared  aperture  components,  i.e., 
the  tilt,  focus,  and  deformable  mirrors.  Thermal  distortions 
of  the  relay  mirrors  affect  only  the  outgoing  wave. 


2.2  OPTICAL  COMPONENT  MODELS 

2.2.1  Tilt  Mirror  Model 

The  tilt  mirror  model  was  based  on  component  infor- 
mation from  Ref.  3 and  was  structured  to  include  the  following 
effects  : 

• Mirror-flexure  hinge  dynamics 

• Voice  coil  actuator  dynamics 

• Mount  compliance  dynamics 

• Laser  induced  vibration. 

Figure  2.2-1  depicts  the  physical  model  and  includes  a listing 
of  the  state  variables,  inputs,  and  disturbances  used  in  the 
mathematical  model. 

The  governing  differential  equations  for  the  pitch- 
tilt  dynamics  are  given  by  Eqs . 2.2-1  to  2.2-4.  A listing  of 
the  relevant  system  parameters  is  shown  in  Table  2.2-1.  The 
equation  format  and  parameters  for  yaw  tilt  (i.e.,  p ) are 
identical. 


x - Fx  + £0cmd  + heb 


(2.2-1) 
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Figure  2.2-1  Tilt  Mirror  Mechanical  Model  and  Variable 

Definition  for  Mirror  Pitch  Motion 
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TABLE  2.2-1 

TILT  MIRROR  MODEL  PARAMETERS 

T-1756 


SYMBOL 

DEFINITION 

^h 

Surface-flexure  hinge  natural  frequency 

Jin 

Mount  natural  frequency 

I ' 

Mirror  inertia/mount  inertia 

T1 

Voice  coil  time  constant 

R 

Mirror  Surface  Radius 

i ! 

L 
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-a. 


( 


T 

£ - 

(0,0,0,0,1/T1) 

(2.2-3) 

. T 

h ■ 

(0,0,0, 0)^,0) 

(2.2-4) 

Aa3  - 

( Ac1 , Ac2 ,0,0,0) 

(2.2-5) 

A^  * 

29  Rc 
s s 

(2.2-6) 

*c2  ’ *s*s 


(2.2-7) 


Where  9 , to,  and  R are  the  elevation  angle,  azimuth  angle,  and 
s s s 

beam  radius  of  the  mirror  surface. 


2.2.2  Focus  Mirror  Model 

The  focus  mirror  model  represents  the  dynamics  of  a 
Cassegrain  beam  expander  with  a movable  secondary  mirror.  The 
mathematical  model  was  structured  as  the  translational  motion 
analog  of  the  tilt  mirror  model,  i.e,,  the  following  effects 
are  considered: 


Mirror-suspension  spring  dynamics 

Voice  coil  actuator  dynamics 

Mount  compliance  dynamics 

Laser  induced  vibration. 
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Figure  2,2-2  shows  the  mechanical  model  of  the  moveable 
secondary  mirror  and  includes  a listing  of  the  variables 
used  in  the  mathematical  model. 


The  governing  state-space  equations  for  the  focus 
mirror  dynamics  are  given  by  Eqs.  2.2-8  to  2.2-12.  Model 
parameters  are  listed  in  Table  2.2-2. 


x ■ Fx  + gu  + hxb  (2.2-8) 


mirror  state  vector: 


t 

X 


x * MIRROR  SURFACE  D IS PLACEMENT 
s 


x « MOUNT  COMP  f.  I A NCR  DISPLACEMENT 
m 

j » NOHMAt.I7.En  VOICE  coil,  FORCE 


INPUT  : 

n - FOCUS  CONTROL  COMMAND 


DISTURBANCES 

x.  - BASE  MOTION  DERIVED  FROM 
VIBRATION  MOD  El, 


TABLE  2.2-2 

FOCUS  MIRROR  MODEL  PARAMETERS 


T-1757 


SYMBOL 


DEFINITION 


u 

s 

U>_ 

m 

m' 

T 


dmax 


Movable  mass-spring  natural  frequency 
Mount  natural  frequency 
Mirror  inertia/mount  inertia 
Voice  coil  time  constant 
Maximum  mirror  displacement  > | xg  | 
Primary  mirror  focal  length 
Primary  mirror  radius 
Zernilce  focus  coefficient 


0 10 


0 0 


-00 


2 

s 


0 


0 


“s”' 


L o 


0 -(w2m’+oo2)  0 

vs  nr 


0 


0 


-u)2m' 


1 

’t 


£T  - (0,0,0,0,1/r) 
hT  « (0,0,0, u2 , 0 ) 


(2. 


(2. 


2-10) 


2-11) 


(2.2-12) 


2.2.3  Deformable  Mirror  Model 


The  deformable  mirror  model  was  structured  to  include 
the  predominant  effects  of  mirror  structural  resonances  and 
spring  rate  nonlinearities.  Generality  has  been  maintained 
in  the  model  in  that  variations  in  structural  resonant  fre- 
quency, damping,  and  cross-coupling  among  the  optical  modes 
implemented  on  the  deformable  mirror  can  be  conveniently 
effected.  The  deformable  mirror  model  was  structured  according 
to  the  following  assumptions: 

• The  mirror  will  be  used  to  implement  the 
optical  modes  listed  in  Table  2.1-1 

• Dynamics  of  each  optical  mode  are  second 
order  with  linear  coupling  to  the  remaining 
modes 

• Provisions  are  made  for  mirror  nonlinear 
. surface  spring  rate 

• Piezoelectric  actuator  response  is 
instantaneous . 

The  first  two  assumptions  imply  that  the  zonal  dynamics  of  the 
mirror  are  dominated  by  the  actuator  structural  dynamics  and 
are  essentially  invariant  over  the  surface  of  the  mirror. 

Linear  coupling  between  modes  implies  a variation  of  the 
DC  gain  of  .the  zonal  dynamics  with  respect  to  actuator  ap- 
plied force  over  the  surface  of  the  mirror;  i.e.: 

Az(x  ,y  s) 

•— Tsl K(xo,yo)G(s)  (2.2-13) 

£LC 

Here  AZ  and  f are  the  surface  displacement  and  actuator 
force,  respectively,  at  actuator  position  (xQ,y0).  G(s) 
is  a unity  DC  gain  transfer  function  which  is  independent  of 
actuator  position.  K(xQ,yo)  is  an  influence  function  which 
depends  on  actuator  location. 
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Allowance  for  nonlinearity  of  the  mirror  surface 
spring  rate  characteristic  is  based  on  experience  with  the 
United  Technologies  deformable  mirror,  as  described  in 
Ref.  4.  The  restoring  force  of  the  mirror  surface  will  be 
modeled  as  a quadratic  function  of  mirror  position  with  ad- 
justable linear  and  quadratic  term  coefficients.  Approxima- 
tion of  the  piezoelectric  actuator  response  by  a unity  gain 
is  a simplifying  assumption  which  should  be  acceptable  for 
the  anticipated  range  of  dither  frequencies.  Actuator 
hysteresis  effects  have  not  been  included. 

The  mathematical  model  which  embodies  these  assump- 
tions is  as  follows: 

z * F(z)2  + Bu  (2.2-14) 

z * Vector  of  mode  strength  and  time 
derivative  of  each  mode 

u = Vector  of  mode  commands 

For  the  kth  mode : 

■ z2k-l  " z2k  (2.2-15) 

z2k  - [“k<z2k-l>]2  z2k-l  - dkz2k  + ^2J12k  Ckjzj 

V Y 

Primary  Dynamics  Coupling 

+ [‘“k(z2k-l)]2uk 

■i  Excitation 

i [U)k(z2x-1)]  “ wkn  + ^k5  !z2k-l 

' ti 


(2.2-16) 


(2.2-17) 


I t 


i 


22 


2 

id.  * Zero  surface  displacement  natural 
u frequency  squared 
* 

w.  ■ Nonlinearity  constant  (related  to 
quadratic  surface  spring  rate 
characteristic ) 


2.2.4  Laser  Phasefront  Model 

In  the  AMD  and  PC  system  simulations  it  is  assumed 
that  the  HEL  output  is  processed  by  a beam  cleanup  system 
that  produces  a uniform  phasefront.  The  phasefront  distur- 
bance models  used  in  the  beam  cleanup  are  described  in 
Section  6.3. 


2.2.5  Thermal  Distortion  Model 


The  thermal  distortion  model  is  based  on  the  results 
for  cooled  mirrors  given  in  Ref.  5.  Two  sources  of  mirror 
deformation  are  considered: 

• Thermal  expansion  effects  (mirror 
surface  thickness  change) 

• Mirror  surface  sag. 

The  first  is  characterized  as  a mapping  of  the  incident 
intensity  function  and  depends  on  nominal  mirror  dimensions, 
material  properties,  beam  incidence  angle,  and  elapsed  time. 
The  second  effect  is  a function  of  the  average  incident 
intensity,  beam  incidence  angle,  and  mirror  geometric  and 
material  properties. 

The  relationship  for  the  phasefront  distortion  due 
to  thermal  expansion  is  (Ref.  5): 

<$l(x,y)  - ~ a^I(x,y)  cos20i(  l-exp(t/Tthermal ) ) (2.2-18) 


I 


I 


V 


p 
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The  corresponding  equation  for  surface  sag  is: 

2 ( 

;2taa0isln01(l-exp(t/Tthermal)) 


V*-y>  ■ — 2aIave 


1H 


(2.2-19) 


where : 


I(x,y)  - Incident  beam  intensity 
I ■ Average  beam  intensity 

a>  V0 

a * Surface  absorption  coefficient 

r ■ Radius  from  center  of  beam 

d **  Beam  Diameter 

0i  * Beam  incidence  angle 

X » Radiation  wavelength 

t ■ Elapsed  time 

Ttherraal  " Mirror  thermal  time  constant 

£-,,S2  “ Mirror  property  dependent  factors 

determined  from  empirical  relation- 
ships in  Ref.  2. 


Here  an  exponential  rise  approximation  has  been  made 
to  the  more  complicated  time  dependence  for  <5^  presented  in 
Ref.  5;  the  same  time  dependence  factor  was  assumed  for  the 
mirror  sag  response. 

Both  Eq.  2.2-18  and  2.2-19  implicitly  include  the 
effects  of  intensity  spread  due  to  beam  incidence  angle. 

For  beam  incidence  angles  other  than  90°,  the  incident 
intensity  profile  on  a mirror  is  the  geometric  projection 
of  a rotationally  symmetric  Gaussian  function  onto  the 
plant  of  the  mirror  surface.  The  mirror  surface  distortions 
produced  by  this  nonsymmetric  intensity  distribution  result 
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in  optical  path  differences  which  produce  the  Gaussian  and 
focus  phasefront  aberrations  shown  by  Eq.  2.2-18  and  2.2-19, 
respectively. 

* 

The  thermallj  induced  aberrations  of  the  following 
optical  elements  were  assumed  to  add  coherently  (Ref.  6): 

• Deformable  mirror 

• Tilt  mirror 

• Relay  mirrors  (4) 

• Beam  expander  primary  and  secondary  mirrors . 

The  combined  effects  of  thermal  expansion  and  mirror 
sag  can  be  decomposed,  into  a focus  mode  aberration  plus  a 
residual  rms  phasefrent  term.  These  two  components  are  sepa- 
rately summed  over  the  set  of  optical  elements  as: 


8 


~ l 

i*l 

<Aal)S! 

U-eXp(-t/Ttherfflali>> 

(2.2-20) 

8 

u$  - l 

rms  i‘;1 

<AVss 

{l-exp(-T/Tthermali)) 

(2.2-21) 

where : 

i,  i_ 

* Focus  mode  strength  of  l component 
due  to  thermal  distortion 

A 1- 

* Residual  phasefront  error  of  i 
component  due  to  thermal  distortion 

ss  * Subscript  indicating  steady  state 
2.2.6  Laser-Induced  Vibration  Model 


A representative  frequency  spectrum  for  laser-induced 
vibration  was  taken  from  Ref.  7.  As  Fig.  2.2-3  shows,  most 
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Laser  Induced  Vibration  Spectrum 
and  Model  Approximation 


of  the  spectral  energy  is  located  in  a band  between  14  and 
50  Hz.  This  spectral  shape  suggests  approximation  by  the 
output  of  an  underdamped  second  order  filter  driven  by 
white  noise.  The  input  noise  variance  and  zero  frequency 
gain  of  the  filter  are  both  set  equal  to  unity;  the  vibration 
"input  gains"  of  the  affected  optical  components  are  adjusted 
to  achieve  a realistic  level  of  component  jitter. 

2.2.7  Relay  Mirror  Jitter  Model 

Figure  2.2-4  shows  the  assumed  relay  mirror  arrange- 
ment; Eqs.  2.2-22  to  2.2-24  define  the  related  mathematical 
model.  The  pitch  and  yaw  plane  dynamics  of  each  mirror  are 
modeled  as  second  order  systems  driven  by  the  output  of  the 
laser  vibration  model.  Due  to  the  spatial  and  temporal  dis- 
tribution of  the  mounting  structure  response,  the  amplitudes 
and  phases  of  the  individual  mirror  responses  will  not  be 
mounting  structure  dynamics,  each  mirror  will  be  assigned  a 
randomly  chosen  vibration  input  gain  (k^)  and  a random  time 
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Figure  2.2-4  Assumed  Relay  Mirror  Layout 


shift  of  the  laser  vibration  model  output.  Geometric  analysis 
shows  that  a mirror  displacement  about  the  pitch  axis  (as 
indicated  in  the  figure)  results  in  a pitch  tilt  aberration  of 
twice  the  displacement  magnitude;  dispalcements  about  the 
yaw  axis  generate  an  amount  of  yaw  tilt  equal  to  the  dis- 
placement . 


Pitch  tilt : 

Yaw  tilt : 

Total  Dynamics : 


0TOTAL  " 2iI10i 
4 

'‘'TOTAL  * i!/1 

si  + 2!i“pt  ®i  + “p^l 


(2.2-22) 

(2.2-23) 

ktVb(Ti}  (2<2"24> 


■ Vibration  "input  gain"  of  ith  mirror 

v^  ■ Vibration  model  output 

t.  * Random  bias  (simulates  re- 
lative phase  of  mirror 
responses) 
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2.2.8  Sensor  Models 


Receiver  Model  (AMD  only)  - The  receiver  is  repre- 
sented by  a fixed  gain.  The  measurement  of  target  return  is 
corrupted  by  internal  noise  of  the  photomultiplier  and  by 
ambient  radiation.  Equation  2.2-25  is  used  to  compute  the 
measured  return. 


m 


(I. 


VKpmt 


N 


pmt 


(2.2-25) 


where : 


I * Measured  return 
m 

- Speckle  corrupted  return 

N » Ambient  radiation  level 
a 

Kpmt  * Photomultiplier  gain 

Npmt  * pb°tomultiPlier  internal  noise 


Shearing  Interferometer  - Reference  8 provides  an 
optical  analysis  and  a thorough  description  of  the  operation 
of  the  shearing  interferometer.  The  net  result  of  the 
optical  analysis  of  Ref.  8 is  Eq . 2.2-26: 


z * He  + n 


(2.2-26) 


z is  a 12  element  vector  of  the  electrical  phases  measured 
by  the  subapertures  of  the  shearing  interferometer.  The  columns 
of  the  H matrix  represent  the  observation  of  a particular  aber- 
ration mode  Cj  ; i.e.,Hij  is  the  electrical  phase  measured  by 
subaperture  i due  to  a unit  value  of  aberration  mode  j.  Off- 
line computation  of  is  described  in  Ref.  8.  c is  the  vector 
of  aberration  mode  strengths  at  the  shearing  interferometer 
due  to  atmospheric  and  thermal  distortion  effects: 
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A 

C = CS  + CA  (2.2-27) 

where  CS,  and  CA  are  the  shared  aperture  optical  component  and 
atmospheric  aberration  modes,  respectively  (see  Fig.  2.1-2). 
n is  a vector  of  zero  mean  white  noise  processes;  the  variance 
of  these  measurement  noise  processes  is  inversely  proportional 
to  the  signal-to-noise  ratio  at  the  shearing  interferometer. 
The  determination  of  this  signal  to  noise  ratio  is  described 
in  Section  2.3.3. 


2.3  ATMOSPHERIC  PROPAGATION  EFFECTS 

2.3.1  Atmospheric  Turbulence  Models 

Frequency  spectra  corresponding  to  the  controlled 
optical  mode  shapes  defined  in  Table  2.1-1  are  generated 
using  an  AFWL-supplied  program,  Ref.  1.  The  spectrum  for 
each  mode  is  dynamically  modeled  as  the  output  of  a linear 
shaping  filter  driven  by  white  noise.  The  computed  turbu- 
lence spectrum  and  filter  spectrum  approximation  for  a tilt 
mode  are  shown  in  Fig.  2.3-1;  Fig.  2.3-2  is  the  shaping  filter 

, realization.  The  output  of  the  turbulence  model  is  also 

delayed  one  transit-time-to-target  (TTTT)  before  it  is  used 

i 

in  the  target  intensity  computation. 

* 2.3.2  Intensity  at  Target 

' The  intensity  distribution  on  the  target  plane  is 

’ , the  squared  magnitude  of  the  Fraunhofer  transform  of  the  ampli- 

I tude  and  phase  distribution  at  the  exit  aperture.  The  ampli- 

tude distribution  is  assumed  to  be  uniform  over  the  aperture 
/ and  the  intensity  distribution  is  computed  as  a function  of 

* f 

, !l  the  phase  distribution  characterized  by  Zernike  coefficients. 
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Figure  2.3-1  Turbulence  Spectral  Representation 
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Figure  2.3-2  Turbulence  Dynamic  Model 
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For  computational  simplicity,  the  exact  distribution  is 
approximated  by  a Gaussian  function  of  the  Zernike  coef- 
ficients, g(c),  as  given  below. 


g(c) 


2tto*L 


1 2 
exp  - 3<1cg3) 


:]eXP[^”  ®XP['t(k 


C3}  j 


(2.3-1) 


where : 
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(2.3-2) 

(2.3-3) 

(2.3-4) 

(2.3-5) 


Pt»  Transmitted  power 
X ■ Laser  wavelength 
Z = Distance  to  the  target 
R » Radius  of  the  aperture 


2.3.3  Target  Effects 


In  both  the  phase  conjugate  and  active  multidither 
systems,  interactions  of  the  incident  laser  beam  with  the 
target  surface  affect  the  optical  feedback  signal.  In  the 
phase  conjugate  model,  the  heating  of  the  target  surface 
affects  its  radiation  level;  for  active  multidither,  target 
surface  shape,  roughness,  and  reflectivity  variations  corrupt 
the  return  wave  with  speckle.  Mathematical  models  of  target 
effects  are  described  in  this  section. 


Speckle  Model  (AMD  only)  - Speckle  is  represented 
as  a normalized  (unity  mean)  multiplicative  noise  which 
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modulates  the  target  return  intensity..  In  general,  the  spec- 
tral content  of  the  speckle  modulation  is  a complicated  func- 
tion of  target  surface  properties  and  target  motion.  Two 
limiting  cases,  however,  provide  useful  guidelines  for  the 
construction  of  a simple  model.  First,  if  the  target  surface 
is  perfectly  specular,  the  return  intensity  is  unmodulated; 
hence,  speckle  could  be  regarded  as  a Gaussian  random  process 
of  unity  mean  and  zero  variance.  The  other  extreme  is  repre- 
sented by  a totally  diffuse  target.  In  this  case,  the  net 
effect  of  speckle  is  modulation  of  the  target  return  by  a 

random  process  of  unity  mean  and  unity  standard  deviation. 

/ 

A reasonable  structure  for  a dynamic  model,  which  would  re-  ■ 
present  a case  between  these  two  extremes,  is  a unity  mean 
Gauss-Markov  process  having  a standard  deviation  proportional 
to  the  diffusivity  of  the  target  surface. 

The  dynamic  system  realization  of  the  simplified 
model  structure  is  a unity-DC-gain  shaping  filter  driven  by 
white  noise  and  a unity  input  as  shown  in  Fig.  2.3-3.  The 
appropriate  frequency  domain  representation  of  the  shaping 


IDEAL 

TARGET 

RETURN 


W (t) 

| 1 

r 

SPECKLE 

SPECKLE  NOISE 

TARGET  RETURN 

1 w 

GW 

A 

Figure  2.3-3  Speckle  Dynamic  Model 


filter  is  determined  from  first  order  filter  cascade  approxi- 
mations to  speckle  spectra  supplied  by  General  Research  Cor- 
poration, Ref.  9.  A cascade  of  three  first  order  filters, 
each  having  a break  frequency  of  2300  Hz,  yields  a good  ap- 
proximation to  that  data  (see  Fig.  2.3-4).  The  input  variance 
is  adjusted  to  yield  a standard  deviation  of  0.3. 


Figure  2.3-4  Speckle  Filter  Approximation 


Target  Radiation  (PC  only)  - The  passive  radiation 
of  a target  is  a complicated  function  of  target  shape,  tem- 
perature (including  heating  effects  of  the  incident  laser 
beam),  surface  characteristics,  and  of  the  ambient  radiation 
level.  The  mathematical  analysis  of  target  passive  radiation 
is  presented  in  Ref.  8;  the  net  effect  of  all  of  the  contribut- 
ing influences  is  the  signal  level  at  the  shearing  interfero- 
meter. For  the  purposes  of  this  study  the  target  effects  model 
will  be  for  a rectangular  target  having  a uniform  radiation  level; 
the  shape  and  radiation  level  determine  the  signal  to  noise 
ratio  of  the  shearing  interferometer  measurements. 
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3. 


MODE  ESTIMATION  IN  ADAPTIVE  OPTICS 


3.1  ACTIVE  MULTIDITHER 

Inaccuracies  in  the  mode  estimates  for  the  AMD  scheme 
are  caused-  by  measurement  noise  corrputing  the  detection  of 
the  reflected  intensity  from  the  target,  as  well  as  by  the 
speckle  effects  due  to  reflections  from  a diffused  target 
surface.  The  mode  estimates  can  be  improved  by  incorporating 
in  the  estimation  algorithm  statistical  information  available 
on  the  turbulence  induced  phase  fluctuations,  measurement  noise, 
and  time  varying  speckle  patterns.  Conventional  (analog  or 
digital)  mode  estimation  algorithms  do  not  use  this  apriori 
statistical  information  explicitly.  A Kalman  filter  uses  a 
priori  statistical  properties  to  arrive  at  the  best  estimate 
of  the  aberrated  phasefront.  In  addition,  the  Kalman  filter 
algorithm  provides  estimates  of  the  full  states  of  the  dis- 
turbances. As  described  in  the  next  section,  optimal  control 
of  the  adaptive  optics  components  utilizes  the  full  state  of 
the  disturbance  process  it  attempts  to  counteract. 

The  received  intensity  is  a function  of  the  phase- 
front  introduced  by  the  turbulence  existing  two  transit  times 
in  the  past;  one  transit  time  each  for  propagation  to  and 
from  the  target.  Processing  (conventional  or  Kalman)  of  the 
received  intensity  only  provides  a two-transit -times  delay 
estimate  of  the  phasefront.  Because  of  signil’ leant  time 
delays  encountered  in  long  range  applications,  such  estimates 
differs  substantially  from  the  present  aberrat Lon , i.e.,  that 
to  be  compensated  by  the  active  optics  elements.  However,  as 
the  Kalman  filter  yields  delayed  estimates  of  the  phasefront 
as  well  as  its  time  derivatives,  the  output  o':  the  Kalman 
filter  along  with  the  statistical  information  on  the  phase 


fluctuations  provide  a basis  for  predicting  the  present  aber- 
ration from  the  delayed  estimates. 

In  this  report,  extended  (because  of  nonlinear  mea- 
surements) Kalman  filtering  algorithms  are  developed  to  esti- 
mate the  Zernike  modes  of  the  phasefront  aberrated  by  atmospheric: 
turbulence.  First,  filtering  and  prediction  algorithms  for 
an  ideal  point  target  are  described.  The  modification  to  a 
Kalman  filter  made  necessary  by  the  time  varying  speckle  pat- 
tern in  the  received  intensity  is  discussed  later. 

The  transfer  functions  derived  for  mode  fluctuations 
discussed  in  Section  2.3.1,  and  the  Gaussian  intensity  distri- 
bution at  the  target  as  a function  of  these  modes  (described 
in  Section  2.3.2  form  the  basis  of  the  Kalman  filter  formulation. 
From  the  pole/zero  values  representing  the  spectrum  of  each  mode 
fluctuation  due  to  atmospheric  turbulence  (Figs.  2.3-1  and  2.3-2), 
a continuous  time  state  variable  equation  is  formulated  as, 

CA  - F CA  + GW  (3.1-1) 

The  state  transition  matrix,  $,  is  then  computed  as 

<f(t)  - eFt  (3.1-2) 

For  a sampling  period,  T,  the  discrete  time  Markov  model  for 
the  mode  fluctuations  can  now  be  written  as 

CAk.  * «(T)CAk--1  + Wk-1  (3.1-3) 

where,  the  covariance,  Q,  of  the  discrete  sequence,  W^,  is 
given  by 
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Q * $(T-t)GQ  ,GT$T(t-T)dt 


0 


where , 


EIGw(t)  wT(t)GT]  - Q 5( t-T ) 


(3-1-4) 


(3.1-5) 


The  mode  strengths  in  the  exit  pupil  phasefront  being 
propagated  to  the  target  is  the  sum  of  the  turbulence  induced 
modes,  CA,  and  the  modes  out  of  the  adaptive  optics  components, 
CC . The  measured  intensity,  I,  is  the  Gaussian  function  (Eq. 
2.3-1)  of  the  exit  pupil  mode  strengths  two  transit  times  in 
the  past,  corrupted  with  zero  mean,  white,  Gaussian  measurement 
noise,  V.  That  is, 


I,  - g( CA,  . , + CC,  , , ) + V,  (3.1-6) 

k — k-delay  — k-delay 1 k 

where  the  delay  equals  twice  the  transit  time.  The  same 
expression  is  being  used  for  received  intensity  or  for  the 
intensity  distribution  at  the  target,  since  multiplying 
factors,  e.g.,  target  reflectivity,  solid  angle  of  collection, 
detector  responsivity , etc.,  can  be  accounted  for  by  choosing 
an  appropriately  scaled  variance  for  the  measurement  noise, 

V^.  Equations  3.1-3  and  3.1-5  constitute  the  "truth"  model 
portion  in  the  evaluation  of  estimation  algorithms  for  turbu- 
lence induced  mode  fluctuations  in  an  AMD  system. 


In  the  Kalman  filter  formulation,  the  mode  estimates 
are  extrapolated  between  the  sampling  intervals  by 

CAk(-)  - <HT)CAk_1  (3.1-7  ) 

The  residual  error  between  the  measured  intensity  and  the  ex- 
pected intensity  is  given  by 

vk  * 1k  “ g(2ik-delay(‘)  + -k-delay 5 (3.1-8) 
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The  estimates  are  then  updated  by 


— k-delay(  + ) 


--k-delay(-) 


+ Vk 


(3.1-9) 


where  Kk  is  the  Kalman  gain  vector  defined  later.  The  esti- 
mates of  the  present  aberration,  which  is  to  be  compensated 
by  the  active  optics  component,  is  then  predicted  by, 


CAk  = $ ( delay ) CAk_deiay(+)  (3.1-10) 

where  <t(delay)  is  the  state  transition  matrix  corresponding 
to  the  turbulence  statistics  over  the  length  of  the  delay, 
i.e.,  Eq.  3.1-2  recalculated  with  the  argument  equalling 
twice  the  transit  time.  Comparing  Eqs.  3.1-6  and  3.1-9,  it 
is  observed  that  the  optimal  prediction  is  simply  optimal 
filtering  in  the  absence  of  measurements. 

Computation  of  the  Kalman  gains  requires  a measure  of 
estimation  error  covariance, 

Pk  « E{(CAk  - CAk)(CAk  - CAk)T|  (3.1-11) 

The  covariance  is  extrapolated  between  the  intervals  by 

Pk(-)  - *<T)Pk_1(+)  *T(T)  + Q (3.1-12) 


where  <f(T)  and  Q are  given  by  Eq . 3.1-2  and  Eq . 3.1-4,  respec- 
tively. The  Kalman  gains  are  computed  as 


Kk+delay 


v-kkv-k + "j'1 


(3.1-13) 


where  R is  the  covariance  of  the  measurement  noise,  V,  and 
the  elements  of  Hk  are  given  by 


IS. 


3CA 


k,i 


CA  - CAk(-) 


(3.1-14) 
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The  error  covariance  matrix  is  updated  by 


V + > - [r  - Wlay  Hk]V-> 

The  necessity  of  dithering  the  modes  for  estimating 
them  is  well  understood  in  a conventional  AMD  scheme.  The 
dither  waveforms  provide  carriers  for  the  fluctuating  mode 
strengths  which  are  recovered  by  demodulation.  It  is  of 
interest  to  the  system  engineer  to  reformulate  the  problem 

\ 

and  substantiate  the  requirement  of  dithers  for  observability 
of  the  states  from  the  measurements.  Computing  the  observation 
matrix,  H^,  for  constant  mode  strengths,  it  can  be  shown  that 
the  pair  (4>,  H^)  results  in  an  unobservable  system  without  dither 
signals.  This  is  to  be  expected,  as  only  the  squares  of  the 
states  enter  the  measurement  equation.  Introducing  dither 
signals  cause  the  observation  matrix  to  be  time  varying  even 
for  constant  state  values,  and  thus  provide  observability. 

From  an  engineering  standpoint,  it  is  desirable  to  keep 
the  dither  frequency  as  low  as  possible  while  providing  a suf- 
ficient data  rate  for  compensation  of  aberration.  In  the  digital 
implementation,  new  information  enters  the  system  for  reestima- 
tion of  the  mode  strengths  only  when  a dither  signal  changes 
sign.  As  all  the  modes  are  coupled  through  the  observation 
equation,  dithering  of  a single  mode  provides  information  on 
all  the  modes.  A time  multiplexed  dither  sequence  is  chosen 
to  provide  a high  data  rate  while  maintaining  a low  dither 
frequency  for  individual  modes,  i.e.,  the  modes  are  dithered, 
one  at  a time,  in  a specified  sequence  while  all  the  modes  are 
estimated  when  any  of  the  modes  is  dithered. 

The  effect  of  time  varying  speckle  pattern  in  the 
received  intensity  is  discussed  in  Section  2.3.3.  As  in  the 
case  of  turbulence  modeling,  a discrete-time  state  variable  model 


i 


t- 

i< 


r 


is  derived  to  correspond  to  the  speckle  spectrum  of  Fig.  2.3-4 
In  the  truth  model  of  the  simulation,  the  point  target  return 
in  Eq.  3.1-5  is  multiplied  by  the  output  of  the  speckle  model, 
qk,  so  that  the  received  intensity  becomes 

h * (1  + V^k-delay  * C£k-delay>  + Vk  t3-1'161 

A filter  design  can  be  constructed  which  estimates 
the  output  of  the  speckle  model  and  thereby  recovers  the  ideal 
point  return  from  the  received  intensity.  However,  any  such 
design  will  assume  a knowledge  of  the  true  speckle  model,  in 
addition  to  requiring  a very  high  data  rate.  For  a target  depen- 
dent phenomenon,  availability  of  such  specific  information  is 
questionable.  In  our  modification  to  the  extended  Kalman  filter 
to  account  for  the  speckle  effects,  only  a knowledge  of  the  rms 
value  of  the  temporal  fluctuations  caused  by  speckle  (and  that 
it  is  a broadband  effect)  is  assumed. 

To  accomplish  the  aforementioned  modification,  Eq.  3.1-15 
is  first  rewritten  as 

Ik  * g^~k-delay  + — k-delay^ 

- ^k  k-delay  + ^k-delay>  + Vk  (3-1'17) 

Note  that  the  first  term  is  the  ideal  point  return,  and  that 
the  second  term  is  broadband.  In  the  modified  extended  Kalman 
filter,  the  second  term  is  lumped  with  the  detector  noise,  V, 
and  the  combined  variance,  Rt , is  given  by 

Rt  * °q  *<£>lc  - CD  * B (3.1-18) 

where  OD  is  the  vector  of  dither  amplitudes,  and  it  is  impli- 
citly assumed  that  the  system  is  performing  near  its  optimum, 
i.e.,  the  deflection  from  the  aimpoint  is  solely  due  to  dither. 
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The  Kalman  gains  in  the  modified  extended  Kalman  filter  are 
computed  with  Rt  replacing  R in  Eq . 3.1-12.  By  appropriately 
increasing  the  effective  measurement  noise  variance  in  the 
computation  of  Kalman  gains,  the  effect  of  speckle  is  diluted 
and,  consequently,  a reduced  emphasis  is  placed  on  the  measure- 
ments . 


3.2  PHASE  CONJUGATE 

In  the  PC  scheme,  inaccuracies  in  the  mode  estimates 
are  caused  by  measurement  noise  corrupting  the  sensing  of  the 
phasefront  of  the  passive  target  radiation.  As  in  the  AMD 
scheme,  the  mode  estimates  can  be  improved  by  incorporating  in 
the  estimation  algorithm,  statistical  information  available  on 
the  turbulence  induced  phase  fluctuations  and  the  measurement 
noise.  The  least  square  fit  algorithm  described  in  Section  3.2.1 
does  not  use  such  apriori  statistical  information;  a least  square 
fit  followed  by  a low  pass  filter  (discussed  in  Section  3.2.2) 
makes  a partial  use  of  the  available  information.  However, 
the  Kalman  filtering  algorithm  uses  all  available  inputs  in  an 
optimal  manner  to  afrive  at  the  best  estimate  of  the  aberrated 
phasefront.  Also,  the  Kalman  filter  provides  the  full  state 
vector  of  the  phasefront  aberrations  required  for  the  optimal  a.1 
control  of  the  active  optics  components.  It  should  be  noted  that 
the  phase  measurement  in  the  PC  scheme  correspond  to  the  present 
aberration  (because  it  corresponds  to  the  return  wavefront  through 
the  turbulence  at  or  near  the  aperture),  which  is  unlike  the  in- 
tensity measurement  in  the  AMD  scheme  which  corresponds  to  the 
delayed  aberration. 

The  discrete  time  Gauss-Markov  model  to  simulate  the 
turbulence  induced  mode  fluctuations  is  given  by  Eq.  3.1-3, 
repeated  below 

CAk  - «(T)CAk-1.+  Ik_1  (3.2-1) 


The  mode  strengths,  CE,  in  the  exit  pupil  phasefront  propagated 
to  the  target  are  the  sum  of  the  turbulence  induced  modes,  CA, 
the  modes  from  the  elements  shared  by  the  outgoing  and  incoming 
path,  CS,  and  modes  not  sensed  by  the  shearing  interferometer, 
CL.  The  modes  sensed  by  the  shearing  interferometer  are  the  sum 
of  CA  and  CS,  and  as  detailed  in  Ref.  8 the  electrical  phase 
measurements,  Z,  out  of  the  shearing  interferometer  are  given 
by 


Zk  - H( CS  + CA)k  + Vk  (3.2-2) 

Equations  3.2-1  and  3.2-2  constitute  the  truth  model  portion 
in  the  evaluation  of  estimation  algorithms  for  turbulence  in- 
duced mode  fluctuations  in  a PC  system. 

3.2.1  Least  Square  Fit 

In  this  method,  the  dynamics  of  the  modes  to  be  esti- 
mated and  the  spectral  properties  of  the  measurement  noise  are 
ignored.  For  least  square  estimation  error,  the  estimate  of 
the  net  aberration  is  given  by 

(6s“T~CA)k  » (HTH)"1HTZk  (3.2-3) 

As  the  name  implies,  this  estimate  of  the  net  aberration  achieve 
the  least  square  fit  to  the  observed  data. 

3.2.2  Least  Square  Fit  Followed  By  Low  Pass  Filter 

In  this  method,  the  exact  dynamics  of  the  modes  to 
be  estimated  are  not  considered;  however,  a spectral  separation 
of  the  modes  and  the  process  noise  is  attempted.  The  mode 
estimates  are  given  by: 
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a(^I+£l)k_1  + (l-a)Uk 


(3.2-4) 


(CS  + CA) 


where  Uk  is  the  least  square  fit  estimate  provided  by  Eq . 
3.2-3.  The  value  of  'a'  is  between  zero  and  unity,  and  is 
determined  from  the  desired  bandwidth  of  the  low  pass  filter. 


3.2.3  Kalman  Filter 
— 7 

The  Kalman  filtering  algorithm  optimally  uses  both 
the  modal  dynamics  and  the  measurement  noise  statistics.  The 
filtering  equations  are  essentially  a repetition  of  Eq . 3.1-6 
through  Eq.  3.1-8,  written  for  the  AMD  scheme.  The  mode  esti- 
mates are  extrapolated  between  the  sampling  intervals  by 


CAk(-)  - <5  (T  )CAjc_1  (3.2-5) 

The  residual  errors  between  the  measured  phases  and  the  expected 
phases  is  given  by 


vk  - Zk  - H ^CA k ( - ) + CSk]  (3.2-6) 

The  estimates  are  then  updated  by 

CAk(+)  - CAk(-)  + Kkvk  (3.2-7) 

where  Kk  is  the  Kalman  gain  vector.  For  computing  the  Kalman 
gains,  the  covariance  of  the  estimation  error  (Eq.  3.1-10) 
is  extrapolated  between  the  intervals  by 

Pk(-)  » <XT)Pk-1(  + ) $T(T)  + Q (3.2-8) 

where  $(T)  and  Q are  given  by  Eq.  3.1-2  and  Eq.  3.1-4,  res- 
pectively. The  Kalman  gains  are  computed  as 

Kk  - Pk(-)HT[^HPk(-)HT  + r]-1  (3.2-9) 
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where,  R is  the  covariance  of  the  measurement  noise,  V.  The 
error  covariance  matrix  is  updated  by 

Pk(x)  * [i  - KkH]pk(-)  (3.2-10) 

In  an  actual  implementation  of  the  Kalman  filter,  the 
covariance,  R,  of  the  measurement  noise,  V,  is  not  accurately 
known  in  computing  the  Kalman  gains  (Eq.  3.2-9).  The  heating 
of  the  target  after  the  laser  is  turned  on  increases  the  signal- 
to-noise  ratio  at  the  shearing  interferometer  output.  An  al- 
gorithm is  described  to  change  the  value  of  R to  adapt  for  a 
change  in  the  signal-to-noise  ratio. 

The  residual  error,  v_k,  after  the  extrapolation 
(Eq.  3.2-5)  is  a Gaussian  random  variable  with  covariance, 

Sk>  given  by 

Sk  » HPk(-)HT  + Rt  (3.2-11) 

where  Rt  is  the  true  covariance  of  the  measurement  noise.  A 
scalar  variable,  b,  is  defined  as 

bk  - vk  HPk(-)HT  + R vk  (3.2-12) 

where  R is  the  assumed  covariance  of  the  measurement  noise  used 
in  the  Kalman  gain  computations.  The  expected  value  of  bk  pro- 
vides a measure  of  the  mismatch  between  the  true  and  the  assumed 
values  of  the  measurement  noise: 


E(V 

H A 

for 

R > 

R » 

Rt 

Rt 

(3.2-13) 

> 1 

R < 

K 

An  adaptive  filtering  algorithm  to  change  the  value  of  R to 
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be  used  in  the  Kalman  gain  computations  is  therefore  formulated 
as 


- Rk  - AR  if  bk  < ^ 

Vl  = Rk  if  ^ \ ^ (3.2-14) 

= Rk  + AR  if  > £ 2 

3.3  EVALUATION  OF  THE  ESTIMATION  ALGORITHMS 

3.3.1  Active  Multidither 

The  Kalman  filtering  algorithms  for  the  AMD  system 
is  evaluated  for  two  ranges,  100  Km  and  500  Km.  To  isolate  the 
effects  of  the  transit  delay,  the  same  signal-to-noise  ratio  of 
20  is  assumed  for  both  the  ranges.  The  aperture  radius  is  as- 
sumed to  be  1 m,  and  the  turbulence  parameters  are  listed  in 
Table  3.3-1. 

To  evaluate  the  effects  of  the  data  rate,  the  system 
is  simulated  at  two  dither  frequencies,  250  Hz  and  500  Hz. 

The  filtering  algorithm  for  an  ideal  point  target  is  described 
by  Eq.  3.1-6  to  Eq . 3.1-14;  the  results  of  the  simulation  are 
listed  in  Table  3.3-2.  The  Strehl  ratio  without  the  dither 
loss  is  intensity  computed  with  estimation  error  only.  In 
all  cases,  very  significant  improvements  in  the  Strehl  ratio 
are  achieved  over  that  with  no  estimation  and  control. 

The  modifications  in  the  filtering  algorithm  due  to 
the  speckle  effects  are  described  in  Eq . 3.1-15  to  Eq.  3.:'  -17. 
The  results  of  the  simulation  for  a speckle  standard  deviation 
of  0.3  and  a triple  corner  frequency  of  2.3  kHz,  are  listed  in 
Table  3.3-3.  Although  there  is  a degradation  in  the  Strehl 
ratio  compared  to  the  corresponding  point  target  case,  there 
is  still  a substantial  imrpovement  over  that  with  no  esti- 
mation and  control. 
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TABLE  3.3-1  ' 

TURBULENCE  PARAMETERS  (AMD) 


T-1904 


ZERNIKE  MODES 

FIRST  CORNER 
FREQUENCY  (Hz) 

ROOT  MEAN 
SQUARE  (um) 

1.  r cos  8,  TILT 

2.25 

1.2 

2.  r sin  9,  TILT 

2.25 

1.2 

3.  2(r2  - 1).  DEFOCUS 

4.5 

0.23 

4.  r2sin  28,  ASTIGMATISM 

2.0 

0.33 

5.  r2cos  20,  ASTIGMATISM 

3.5 

0.33 

TABLE  3.3-2 


AMD  WITH  POINT  TARGET 


T-1905 


TABLE  3.3-3 

AMD  WITH  SPECKLE  RETURN 


STREHL  RATIO 


T-190 


WITHOUT  WITH  ESTIMATION  AND  IDEAL  CONTROL 

EST I MAT  I ON  WITH  DITHER  LOSS  WITHOUT  DITHER  LOSS 
C ’ AND  C0KTR0L  " 250~  Hz  ’ 500  ~Hz“' 1 250  Hz'  500  Hz 


100 

.315 

.668 

.690 

.814 

.838 

500 

. 315 

.600 

.550 

1 

.689 

.636 

3.3.2  Phase  Conjugate 

The  filtering  algorithms  for  the  PC  system  are  evaluated 
for  a range  of  5 Km.  The  aperture  radius  is  assumed  to  be 
0.25  m;  the  turbulence  parameters  are  listed  in  Table  3.3-4. 

TABLE  3.3-4 

TURBULENCE  PARAMETERS  (PC) 


ZERNIKE  MODES 


FIRST  CORNER 
FREQUENCY  (Hz) 


fl 

r cos  0 , 

TILT 

2 . 

r sin  0, 

TILT 

3. 

2( r2-l ) , 

DEFOCUS 

4 . 

r2sin  20  , 

ASTIGMATISM 

5. 

r2coe  29 , 

ASTIGMATISM 

7.5 


.5 


15.0 


11.6 


ROOT  MEAN 
SQUARE  (um) 


0.21 


0.3 


0.3 


The  shearing  distances  are  0.017  m in  the  X direction, 
and  0.057  m in  the  Y direction.  The  system  is  simulated  for 
three  rms  levels  of  electrical  phase  measurement  noise  at  the 
output  of  the  shearing  interferometer,  0.1  rad,  0.2  rad  and 
0.3  rad,  The  results  of  the  simulations  plotted  in  Fig.  3.3-1 
are  all  for  a data  rate  of  1 kHz. 

As  the  LSF  makes  no  attempt  to  filter  the  measurement 
noise,  at  a high  noise  level  the  performance  of  L&F  is  actually 
worse  than  attempting  no  corrections.  If  the  bandwidth  of  the 
LPF  is  too  narrow  relative  to  turbulence  bandwidth,  and  the 
noise  level  is  very  low,  the  performance  of  LSF-LPF  is  seen  to 
be  worse  than  simple  LSF,  because  the  LPF  is  filtering  out  the 
true  signal.  As  expected  the  performance  of  the  Kalman  filter 
is  superior  at  all  noise  levels. 


RMS  ELECTRICAL  PHASE  MEASUREMENT  NOISE  (rad) 

1 1/SNR) 

Figure  3.3-1  Phase  Conjugate 
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4. 


OPTIMAL  CONTROL  IN  ADAPTIVE  OPTICS 


4.1  INTRODUCTION 

The  goal  of  an  adaptive  optics  control  system  is  to 
continuously  null  disturbances  arising  from  such  effects  as 
turbulence,  vibration,  and  thermal  distortion.  Ideally,  an 
active  optics  component  such  as  a focus  system  would  counter- 
act perfectly  the  focus  disturbances  of  the  system  environment. 
Practically,  such  realities  as  finite  response  time  and  struc- 
tural mode  exitation  preclude  the  possibility  of  perfect  dis- 
turbance nulling.  Modern  control  theory  recognizes  the  non- 
ideal characteristics  of  active  components  and  offers  a realistic 
yet  ambitious  approach;  the  goal  of  an  optimal  controller  is 
to  minimize  the  difference  between  a disturbance  state  that 
degrades  system  performance  and  the  active  element  response 
that  attempts  to  counteract  it. 

This  chapter  will  deal  with  the  application  of 
optimal  control,  distinct  from  the  estimation  problem,  in 
adaptive  optics.  Section  4.2  outlines  the  formulation  of 
the  optimal  disturbance  nulling  control  problem;  quadratic 
performance  indices  used  in  the  A/0  control  synthesis  are 
derived  in  Section  4.3.  The  inner  control  loops  used  to 
extend  the  control  bandwidths  of  the  tilt  and  focus  optics 
are  described  in  Section  4.4.  Application  of  optimal  dis- 
turbance nulling  control  to  the  deformable  mirror  required 
an  unusual  approach  which  is  summarized  in  Section  4.5.  Some 
preliminary  results  which  demonstrate  the  advantages  of  the 
modern  control  approach  over  candidate  classical  designs  are 
provided  in  Section  4.6;  additional  control  performance  re- 
sults are  presented  for  the  integrated  estimation-control 
system  in  Section  5. 
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4.2  OPTIMAL  DISTURBANCE  NULLING  CONTROL 

The  general  structure  of  an  optimal  disturbance 
> nulling  control  system  is  shown  in  Fig.  4.2-1.  The  system 

consists  of  a linear  dynamic  model  of  the  disturbance  to  be 
nulled,  a matrix  of  feedforward  gains  which  operates  on  the 
disturbance  states  to  compute  the  control  , and  an  active 
element  having  a feedback  controller.  In  an  adaptive  optics 
system  the  disturbance  model  might  represent  the  dynamics  of 
a turbulence-induced  aberration  while  the  controllable  element 
would  be  an  active  component  such  as  a tile  mirror.  Here,  ideal 
! estimation  has  been  assumed  in  accordance  with  the  Separation 

Theorem  (Section  1.2);  all  of  the  disturbance  model  states  are 
available  to  the  feedforward  controller. 

The  optimal  structure  shown  in  Fig.  4. .2-1  utilizes 
all  available  knowledge  of  the  dynamics  of  both  the  distur- 
bance model  and  the  controllable  element,  as  well  as  all 
current  information  on  the  disturbance  state  in  computing 
the  controls  u1  and  u.  More  specifically,  the  state  vector 


i 

4] 


DISTURBANCE 

MODEL 


FEED 

FORWARD 

CONTROLLER 


CONTROLLABLE 
ELEMENT  WITH 
FEEDBACK  CONTROLLER 


Figure  4.2-1  Disturbance  Nulling  Control  Structure 
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of  the  disturbance  model  provides  the  present  value  of  the 
disturbance  plus  states  related  to  the  time  derivatives  of 
the  disturbance.  CJsing  the  full  state  of  the  disturbance, 
the  control  u.^  can  be  computed  not  only  to  counteract  the 
present  value  of  the  disturbance  but  also  to  anticipate 
(in  part)  the  future  of  the  process  through  th6  derivative 
states.  It  will  be  shown  in  Section  4.6  that  this  anticipa- 
tion effect  can  lead  to  a significant  alleviation  of  sample 
rate  and  controllable  element  bandwidth  requirements. 

The  formulation  and  solution  of  the  mathematical 
problem  to  determine  K and  K1  car  be  facilitated  by  re- 
drawing the  structure  of  Fig.  4.2-1  into  an  optimal  regulator 
format,  as  shown  in  Fig.  4.2-2.  The  augmented  plant  consists 
of  the  dynamics  of  the  controllable  element  and  of  the 

q-30081 

AUGMENTED 


disturbance;  the  augmented  regulator  includes  both  the  feed 
forward  and  feedback  controllers.  The  optimal  regulator 
problem  then  specifies  that  the  control  gains  are: 

K1  - R"Vp  (4.2-1) 


where : 

/>.rp  rp  | rp 

K1  - (KXJK^) 

R"1  » Inverse  of  control  weighting  matrix 

G * Control  effectiveness  of  augmented  system 

A 

P * Solution  of  Ricatti  Equation  for  the 
augmented  system 


Following  Ref.  2,  the  mathematical  description  of 
the  augmented  system  is  given  by  Eqs.  4.2-2  to  4.2-4. 


;T  - C*V> 

•Arn 

Gx  » (G  jO) 

- “f  ! o ~ 

F - J 

0 j A 

with  the  performance  index: 

J - / SuTRu  + (£+£)T  Q(l+£)jdt 

» JQ  ) uTRu  + (HTx+CTz)T  Q(HTx+CTz) jdt 
- / juTRu  + STQ£fdt 


(4.2-2) 

(4.2-3) 

(4.2-4) 

(4.2-5) 

(4.2-6) 

(4.2-7) 


where : 


i 


4 'f 


1 T -1  T 

Q J QH(HXH)  1C1 

L(QH(HaH)  iCir  i C(HaH)  xCl . 

Q ! Ql 


The  time  invariant  Ricatti  equation  for  the  augmented 
system  becomes 


A A Am  A A A T A A 

o » PF  + FaP  - PGR“AG*P  + Q 


(4.2-10) 


Partitioning  P according  to: 

n l m 


IT* 
p i P 

p i *21 

*21 ! 22 


(4.2-11) 


The  Ricatti  equation  can  be  decomposed  into  Eqs . 4.2-12  to  4.2-14 


PF  + FTP  - PGR‘1GTP  + Q 

(4.2-12) 

P21(F  - GR-1GTP ) + atp21  + Q21 

(4.2-13) 

P22A  + ATP22  - P21GR-1GTpT1  + Q22 

(4.2-14) 

Note  that  Eq.  4.2-12  is  the  Ricatti  equation  corresponding  to 
the  optimal  regulation  of  the  controllable  element  by  itself; 
i.e.,  the  disturbance  dynamics  do  not  affect  this  partition 

A 

of  the  P matrix.  Equations  4.2-12  to  4.2-14  can  be  solved 
sequentially;  for  convenience  Eq.  4.2-13  can  be  rewritten: 


P21Fcl  + A P21  ” “ Q21 


(4.2-15) 

where  F ^ is  the  closed  loop  F matrix  of  the  active  element. 
The  optimal  control  gains  are  given  by: 

(4.2-16) 


KT  - - R_1GTP 
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(4.2-17) 


T -ITT 

Ki  - -R  G pL 


Equations  4.2-15  and  4.2-17  indicate  that  the  optimal 
feedforward  controller  is  determined  by  the  closed  loop  dynamics 
of  the  controllable  element  (Fc^),  the  dynamics  of  the  distur- 
bance (A),  the  state  penalty  matrix  (Qg^),  and  the  control 
penalty  matrix  (R).  Optimal  regulator  theory  insures  that 
the  controller  derived  from  the  augumented  system  analysis 
will  yield  the  minimum  cost  J consistent  with  the  weightings 
and  dynamics  involved. 


4.3  PERFORMANCE  INDICES  FOR  ADAPTIVE  OPTICS 
CONTROL  DESIGN 

In  optimal  control  analysis,  the  success  of  a given 
system  design  is  measured  in  terms  of  scalar  "cost  functional" 
or  performance  index.  For  a regulator  or  tracking  system  the 
performance  index  generally  reflects  the  integration  of  a 
weighted  sum  of  the  system  output  error  and  the  corresponding 
control  effort  over  a specified  period  of  time.  Once  the  per- 
formance index  has  been  chosen,  the  optimal  control  of  the 
system  is  determined  by  finding  the  control  time  history  which 
minimizes  this  index. 

This  section  outlines  the  derivation  of  quadratic 
performance  indices  (PI)  for  candidate  high  energy  laser 
adaptive  optics  systems.  Physical  considerations  lead  to 
an  exponential  weighting  function  for  the  phasefront  aber- 
rations of  the  system;  approximation  by  series  expansion 
reduces  the  exponential  function  to  a quadratic  form.  A 
quadratic  form  is  desirable  since  linear  optimal  disturbance 
nulling  synthesis  is  then  possible.  The  guidelines  for  quadra- 
tic PI  construction  vary  significantly  between  the  active  multi- 
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dither  (AMD)  and  the  phase  conjugate/beam  cleanup  concepts; 
hence,  the  discussion  is  divided  between  the  two  system  types. 

4.3.1  Development  of  the  State  Penalty  Matrix,  q 

The  performance  goal  of  all  of  the  candidate  adaptive 
optics  systems  is  clearly  defined  on  physical  grounds:  the 
system  should  maximize  the  intensity  of  incident  radiation 
at  a point  fixed  to  the  target.  The  most  precise  weighting 
function  would  therefore  be  the  difference  between  the  maximum 
achievable  target  intensity  and  that  actually  produced,  i.e. 


f (c1,c2 ....  Cg) 


I 

max 


I ( c^ , c2  < ^3  > • * • Cg ) 


(4.3-1) 


Conversion  of  Eq . 4.3-1  into  a Q matrix  requires  approximation. 
If  the  c^’s  are  all  small  (as  they  would  be  if  successfully 
controlled),  I (c^ . Cg , , . . . c^)  can  be  approximated  bv  a series 
expansion.  The  general  form  of  the  Taylor  series  for  I is: 


^C1  ’ c2  ’ c3  ’ * * 


V 


I(c,  , Cn  i • ■ . C.  ) 

J'0  ^0  J0  50 

1 5 5 ^ 2 y 

+ lft-  + jl, 

(4.3-2 ) 


Ac  . + . 
3 


Here  the  partial  derivatives  are  evaluated  at  the  maximum 
intensity  values  of  the  ci's,  signified  by  ci  . If  terms  up 
to  second  order  are  retained,  the  approximate°intensity  degra- 
dation function  becomes: 


f ( Ac^ , Ac2 , . . . Ac  ) 


32I 


3c j 


AciAc^ 


(4.3-3) 
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Phase  Conjugate /Beam  Cleanup  Q Matrix  - The  maximum 
target  intensity  for  the  phase  conjugate  and  beam  cleanup 
systems  is  achieved  when  all  of  the  c^'s  are  zero.  At  this 
point  all  first  partial  derivatives  of  Eq . 4.3-1  are  zero 
so  that  the  intensity  degradation  function  reduces  directly 
to  the  desired  quadratic  form: 


f(c 


1 ’ 2 


C7> 


iff  d2I 

”2  ii1  3ci3cj 


AciAcj 


(4.3-4) 


which  is  the  scalar  equivalent  of  the  matrix  form: 


f (c^ .Cg  » . . • c^) 


T 

Ac  Q 


pc/bcA- 


(4.3-5) 


with 


qu 


D2! 


5ci9cj 


(4.3-6) 


Again,  Eq . 4.3-5  expresses  the  reduction  of  target  intensity 
as  a function  of  the  deviations  of  the  ci’s  from  null,  i.e., 
the  desired  maximum  intensity. 


i 


Since  the  c.'s  enter  the  target  intensity  function, 
1 2 2 2 

Eq.  2.3-1,  as  squares,  i.e.  I * ICc-pCg,...  Cg),  all  of  the 
second  derivatives  also  vanish 


32I 


3ci3cJ 


0 for  i^j 


(4.3-7) 


Q,  therefore,  is  a diagonal  matrix  with: 


Differentiating  Eq . 
modes  are: 


1 92I 


2.3-1  the  q^'s 


(4.3-8) 

for  the  tilt  and  focus 


!> 
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q22 

q33 


sf  i 

3C? 


32I 
— ? 
3C2 


32 1 
— 2 
3c3 


‘max  4 


'max  4 


‘max  3 


( 4 . 3-9 ) 

(4.3-10) 

(4.3-11) 


The  <li;L's  for  the  remaining  modes  are  proportional  to  deriva- 
tives of  the  Strehl  ratio: 


*44 


q55 


I a 
max  6 


(4.3-12) 


Physically,  the  q^,  q22  ’ q33  elements  weight  intensity 
degradation  due  to  beam  tilt  and  focus  plane  displacements; 
q44  and  q55  weight  the  relative  beam-broadening  effects  of 
the  astigmatism  modes. 

Active  Multidither  Q Matrix  - When  dither  signals  are 
introduced,  the  maximum  achievable  intensity  no  longer  corresponds 
to  the  stationary  point  of  the  Gaussian  intensity  function.  if 
the  phasefront  is  dithered  with  square  waves,  as  it  would  be  in 
a digital  application,  the  desired  intensity  is  given  by  Eq. 

2.3-1  evaluated  at  the  nominal  dither  amplitude,  i.e. 

c . •*  d.  , i°l,  ...  5 i (4.3—13) 

0 

\ 

In  * I(c-  ,c0  , ...  c,  ) (4.3-14) 

0 10  20  d0 

Since  the  nominal  oporating  point  is  not  at  the  peak 
of  the  Gaussian  function  all  of  the  coefficients  of  the  first 
order  term  of  Eq.  4.3-3  are  nonzero.  Neglecting  the  second 
order  terms  of  Eq . 4.3-3,  the  intensity  degradation  function 
becomes : 
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The  "0"  subscript  indicates  the  nominal  magnitude. 
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4.3.2  Development  of  the  Control  Penalty  Matrix,  R 


While  the  Q matrices  were  derived  from  explicit 
system  goals  and  physical  arguments,  the  R matrices  will  be 
constructed  to  reflect  the  estimated  relative  cost  of  control 
effort,  and  to  maintain  performance  limits  within  the  capa- 
bilities of  the  adaptive  optics  components.  The  rationale 
for  development  of  an  R matrix  depends  on  the  particular 
component  involved  and,  in  the  present  case,  will  necessarily 
be  a result  of  the  design  process.  In  Section  4.4  the  R's 
for  the  tilt  and  focus  mirrors  are  determined  in  designing 
inner  control  loops  for  these  components.  The  R matrix  for 
the  deformable  mirror  is  determined  in  a rather  unconventional 
manner;  Section  4.5  summarizes  this  design  approach. 


4.4  . INNER  LOOP  DESIGN  - TILT  AND  FOCUS  MIRRORS 

In  this  section  the  design  of  the  inner  control  loops 
for  the  tilt  and  focus  mirrors  is  described.  Since  the 
tilt  and  focus  mirror  dynamic  models  are  rotational  and  trans- 
lational analogs  of  each  other,  the  detailed  analysis  presented 
here  is  restricted  to  determining  the  feedback  gains  for 
one  of  the  tilt  mirror  axes.  The  feedback  gains  for  the  focus 
mirror  can  be  obtained  by  a procedure  identical  to  that  used 
for  the  tilt  mirror. 

4.4.1  Control  Objectives  - Performance  Indices 

The  control  objectives  for  the  tilt  mirror  are 

• Extend  nir *or  surface  response  bandwidth  to 
1000  Hz  (*<5280  rad/sec) 

• Maintain  a stable  mount  compliance  mode  response. 
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are : 


The  corresponding  specifications  for  the  focus  mirror 


* 

\ 


• Extend  secondary  mirror  position  response 
bandwidth  to  500  Hz 

• Maintain  a stable  mount  compliance  mode  response. 

These  specifications  reflect  a desire  to  improve  hardware 
response  bandwidths  to  current  or  soon-expected  standards , 
as  well  as  a concern  for  potential  instabilities  of  struc- 
tural (mount  compliance)  modes  that  may  result  from  this  in- 
creased bandwidth. 

Since  the  tilt  and  focus  mirrors  are  single  input - 
single  output  devices  (the  two  tilt  axes  respond  independently) 
the  R matrix  of  the  quadratic  performance  index  reduces  to  a 
scalar  which  is  balanced  against  a single  diagonal  element  of 
the  Q matrix  derived  in  Section  4.3.  The  mathematical  results 
will  show  that  the  bandwidth  of  the  closed  loop  system  is 
determined  by  the  ratio  q/r,  q and  r being  the  state  and  control 
weighting  scalars,  respectively.  With  q fixed  by  the  results 
of  Section  4.3,  r is  chosen  to  achieve  the  specified  response 
bandwidth . 

4.4.2  The  Solution  Technique 

Determining  the  optimal  feedback  gains  for  a time  in- 
variant regulator  involves  solution  of  the  steady-state  Ricatti 
equation.  In  the  case  of  the  tilt  mirror,  the  coefficient 
matrices  for  the  Ricatti  equation,  Eq . 4.4-1,  are  given  by 

1 1 

I Eqs.  4.4-2  to  4.4-6. 
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(4.4-6) 


There  are  a number  of  alternative  methods  for  solving 
Eq.  4.4-1  , such  as  integrating  the  Ricatti  differential  equation 
to  steady  state  or  using  various  numerical  iteration  techniques. 
The  technique  used  in  the  present  study  was  the  singular  per- 
turbation method  as  developed  by  Sannuti  and  Kokotovic  (Ref.  10). 
Singular  perturbation  takes  advantage  of  the  characteristic 
frequency  disparity  between  the  mirror-flexure  hinge  dynamics 
and  the  mount  compliance  and  actuator  responses  to  effectively 
reduce  the  order  of  the  dynamic  system.  A complete  discussion 
of  the  technique  would  be  inappropriate  here  although  some  of 
the  major  features  should  be  mentioned.  (The  interested 
reader  is  encouraged  to  examine  Ref.  10.)  The  technique: 


• Reduces  the  order  of  the  Ricatti  equation  down 
to  that  of  the  dominant  slow  dynamics;  the  fast 
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dynamics  (such  as  tiie  mount  compliance  mode 
dynamics)  ape  considered  as  a perturbation 

• Control  gains  are  obtained  as  a two-term 

perturbation  expansion.  The  leading  order 
term  of  the  expansion  controls  the  dominant 
slow  dynamics;  the  second  term  maintains  stabi 
lity  of  the  fast  dynamics  and  "patches  up" 
the  control  of  slow  dynamics  to  account  for 
coupling  of  the  fast  modes. 


In  this  particular  design,  singular  perturbation  reduces 
the  fifth  order  system  to  the  second  order  dynamics  of  the  mirror- 
flexure  hinge  mode  (i.e.,  x1  and  Xg);  the  mount  compliance  and 
actuator  dynamics  are  considered  as  a perturbations.  The  re- 
duced order  Ricatti  equation  is  simple  enough  that  it  can  be 
solved  analytically.  The  equations  for  the  second  term  of  the 
perturbation  expansion  for  the  gains  can  also  be  solved  analyti- 
cally; the  resulting  gains  are  given  by: 


T T T 

K1  * Kq  + eKA  (4.4-7) 

( b , /2b  ,0,0,0)  + e(-  wmt1/2b,uimT1b,0,-:~,a)mT;L/2b) 
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(4.4-10) 


The  first  terra  of  Eq.  4.4-8  feeds  back  position  and  rate  of 
the  mirror  surface  motion,  thereby  extending  the  mirror  band- 
width and  improving  the  damping  of  the  mirror-flexure  hinge 
mode.  The  second  term  is  of  small  magnitude  compared  to  the 
firsf,  since  it  is  multiplied  by  a small  number,  e.  The  first 
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two  elements  of  the  second  term  "patch  up”  the  mirror  position 
and  rate  terms  to  account  for  the  destabilizing  effects  of 
mount  compliance  mode  coupling  and  actuator  lag.  Mirror  posi- 
tion feedback  is  decreased  slightly;  rate  feedback  is  increased 
slightly;  the  net  effect  is  a small  increase  in  the  damping 
of  the  mirror  surface  dynamics.  The  fourth  element  of  the 
second  term  feeds  back  a control  proportional  to  the  rate 
of  the  mount  compliance  mode.  The  metastable  (undamped) 
dynamics  of  this  mode  became  unstable  as  a result  of  the  Kq 
feedbacks;  a small  amount  of  rate  feedback  wilJ  restabilize 
it.  Finally,  the  bandwidth  of  the  actuator  is  increased 
slightly  by  the  fifth  element  of  the  second  terra. 

The  nacural  frequency  and,  hence,  the  bandwidth  of 
the  reduced  order  dynamics  (second  order  mirror-flexure  hinge 
dynamics  and  gains)  is: 

“2  " ■%  * “hb  (4.4-11) 

- go^n/i  +(!)  (4.4-12) 

w - uhvr  +(Yj  (4.4-13) 

* Mirror  response  bandwidth 

Using  Eq.  4.4-13,  a bandwidth  specification  can  be  readily 
converted  to  a q/r  ratio: 

| .(JL)4  . ! (4.4-14) 

Equations  4.4-8  to  4.4-10  can  be  modified  to  yield 
gains  for  the  focus  mirror  by  a few  simple  changes  in  notation: 
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In  summary,  a controller  has  been  designed  which  meets 
the  specifications  listed  in  Section  4.4.1.  The  solution  tech- 
nique used  yielded  the  optimal  regulator  gains  in  closed  form, 
as  well  as  approximate  expressions  for  the  closed  loop  bandwidth 
Using  the  closed  form  expressions,  new  gains  can  be  readily 
computed  to  accommodate  changes  in  mirror  parameters  or  control 
specifications . 


4.5  DEFORMABLE  MIRROR  CONTROLLER  DESIGN 

In  a number  of  practical  design  situations,  synthesis 
of  an  optimal  control  system  is  hampered  by  the  presence  of  an 
inalterable  dynamic  element.  In  flight  vehicle  control  engineer- 
ing, for  example,  the  lag  dynamics  of  '.he  control  surface  actua- 
tors may  be  inalterable  due  to  the  lack  of  a feedback  measurement; 
i.e.  the  actuator  dynamics  states  are  unobservable.  In  man- 
vehicle  control  analysis,  the  dynamics  of  the  neuro-muscular 
lags  of  the  pilot  are  fixed;  one  cannot  close  an  inner  Loop 
around  the  pilot  himself.  In  the  control  of  higher  order  aber- 
rations in  adaptive  optics  a deformable  mi? ror  is  used  which 
has  inalterable  dynamics,  i.e.  no  measurement  of  the  mirror 
figure  is  available  to  construct  an  inner  feedback  loop. 

The  inalterable  dynamics  of  the  deformable  mirror 
call  fcr  a unique  approach  to  the  optimal  disturbance  nulling 
control  problem;  this  approach  and  the  resulting  feedforward 
controller  design  are  described  in  this  section. 


4.5.1  Pseudo-Optimal  Control  Design 


An  approach  to  optimal  control  synthesis  for  systems 
having  inalterable  elements  is  pseudo-optimal  control.  The 
pseudo-optimal  control  concept  is  to  identify  a dynamic 
subsystem  within  a desired  optimal  system  structure  that  can 
be  replaced  by  the  dynamics  of  fne  inalterable  elements.  After 
the  inalterable  element  is  substituted  into  the  role  of  the 
dynamic  subsystem,  the  remainder  of  the  desired  system  control 
is  optimized  about  that  starting  point.  This  design  methodology 
yields  the  "best”  control  system  which  can  exist  while  also  in- 
corporating the  inalterable  element. 

The  application  of  the  pseudo-optimal  design  philos- 
ophy to  control  of  a deformable  mirror  is  illustrated  by  a 

comparison  of  Fig.  4.5-1  with  Fig.  4.5-2.  The  design  goal  is 

T 

to  find  the  best  gain  matrix  to  be  used  with  the  existing 
mirror  dynamics.  To  this  end,  the  dynamics  of  the  deformable 
mirror  are  cast  in  the  role  of  tne  closed  loop  element  encir- 
cled in  Fig.  4.5-2.  As  Fig.  4.5-1  and  4.5-2  suggest,  the 
pseudo-optimal  control  synthesis  consists  of  finding  a fic- 
titious regulator  problem  havxng  the  deformable  mirror  dynamics 
as  the  final  closed  loop  solution,  and  then  solving  the  remainder 
of  the  optimal  disturbance  nulling  control  problem.  The  general 
mathematical  procedure  to  be  used  in  pseudo-optimal  control 
design  for  the  deformable  mirror  is  outlined  in  the  fol- 
lowing section. 

/ 

4.5.2  The  Mathematical  Design  Procedure 

T 

The  computation  of  the  optimal  gains  K1  to  be  used 
with  the  deformable  mirror  is  considered  as  two  subproblems: 

• Construction  of  an  appropriate  fictitious 
regulator  problem 
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Figure  4.5-1  Deformable  Mirror  Controller 


Figure  4.5-2  Optimal  Model  Following  (Nulling) 

Control  Structure 


• Design  of  the  feedforward  controller. 

The  first  subproblem  is  an  involved  process  which 
requires  some  reader  familiarity  with  asymptotic  (perturbation) 
techniques.  A detailed  solution  of  the  fictitious  regulator 
problem  is  presented  in  Ref.  11;  the  major  points  of  the 
design  process  will  be  summarized  here: 


• Mirror  dynamics  are  characterized  as  a set 
of  weakly  coupled  second  order  systems  (see 
Section  2.2,3) 
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• Mirror  dynamics  and  associated  Ricatti 
equations  are  decomposed  through  asymptotics 

- Primary  second  order  dynamics  become 
uncoupled 

- Cross-coupling  effects  are  treated  as  a 
perturbation  to  primary  dynamics 

• Fictitious  regulator  is  developed  for  the 
decomposed  mirror  model. 


The  result  required  from  the  fictitious  regulator 
problem  is  the  R-1  matrix  (inverse  of  the  control  weighting 
matrix).  Using  the  asymptotic  approach,  the  elements  of  the 
R-1  matrix  were  obtained  in  closed  form;  formulae  for  the 
elements  are  given  by  Eq . 4.5^1  and  4.5-2: 


Rnm  “ :r“C2(l  “ C2)  (diagonal  elements)  (4,5-1) 
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where : 

£,oj  are  the  damping  ratio  and  natural  frequency  of 
the  primary  second  order  dynamics 

fj^  ,f,  are  the  position  and  rate  coupling  terms 

xmn  ^mn  of  primary  mode  n into  primary  mode  m 

With  Eq . 4.5-1  and  4.5-2  the  feedforward  controller 
can  be  designed  directly;  the  fictitious  regulator  problem 
need  not  be  solved  for  every  design.  Repeating  results  from 
Section  4.2: 


T -ITT 

Ki  - - R gV21 


(4.5-3) 
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f 

where  P21  is  derived  from  the  mirror  and  disturbance  dynamics 
according  to: 

P21F  + ATp21  *-«21  <4-S'4> 

) 

where  F and  A are  the  dynamics  matrices  of  the  mirror  and 
disturbance,  respectively. 

In  summary,  the  essential  steps  in  designing  a pseudo- 
optimal  controller  or  the  deformable  mirror  are: 

~1 

• Compete  the  elements  of  the  R matrix  of 
the  fictitious  regulator  using  Eq.  4.5-1 
and  4.5-1 

• Determine  P21  from  Eq.  4.5-4  using  the 
mirror  F matrix  and  the  A matrix  of  the 
disturbance  models 

T 

• Calculate  the  feedforward  gains  K 

from  Eq.  4.5-3.  1 

T 

A final  scalar  adjustment  of  the  matrix  is  required; 
the  justification  of  this  adjustment  is  given  in  Section  4.6. 

4.5.3  Description  of  the  Deformable  Mirror  Controller. 

T 

The  controller  has  primary  and  secondary  functions. 
The  primary  function  of  the  controller  is  to  command  a linear 
combination  of  the  states  of  a disturbance  model  to  its  cor- 

1 

responding  primr.ry  dynamics  in  the  deformable  mirror.  The 

• secondary  function  is  to  crossfeed  linear  combinations  of  the 

remaining  disturbance  model  states  to  suppress  cross-coupling. 

. * ' These  two  control  functions  are  depicted  for  a two  disturbance 

i j 

system  in  Fig.  4.5-3  and  4.5-4. 

i 

« 'r 

' !! 
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The  primary  and  secondary  functions  can  be  readily 
extended  to  an  arbitrary  number  of  coupled  second  order 
systems,  i.e.  the  dynamics  of  the  deformable  mirror  as  des- 
cribed in  Section  2.2.3.  This  extention  is  described  in 
Ref.  11. 


4.6  PRELIMINARY  PERFORMANCE  EVALUATION 

Performance  and  sensitivity  studies  of  modern  control 
designs  for  adaptive  optics  components  are  presented  in 
this  section.  Again,  the  estimation  problem  has  been  separated 
from  the  control  design  (disturbance  states  are  available  direct- 
ly from  the  truth  model)  so  that  control  performance  can  be 
isolated  from  all  other  effects  that  may  influence  system 
performance.  Studies  were  conducted  for  the  tilt  channel  and 
the  deformable  mirror  channel;  since  the  focus  channel  has 
directly  analogous  disturbance  and  mirror  dynamics,  perfor- 
mance results  similar  to  those  for  the  tilt  channel  are 
implied. 


The  following  studies  will  be  presented  here: 


• Tilt  Channel  (and  Focus) 

- Disturbance  nulling  performance  compared  to 
that  of  a classical  design 

- Sensitivity  of  optimal  controller  performance 
to  tilt  mirror  bandwidth 

- Sensitivity  of  optimal  controller  performance 

to  truth  model/design  model  mismatches 

• Deformable  Mirror  Channel 

- Disturbance  nulling  performance  including 
mirror  modal  cross-coupling  effects 

- Comparison  of  performance  with  classical  inverse 
filter  controller. 
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4.6.1  Tilt  Channel  Performance/Sensitlvlt 


Figures  4.6-1  and  4.6-2,  respectively,  show  block 
diagrams  of  an  optimal  disturbance  nulling  controller  and  a 
classical  direct  measurement -command  controller  for  the  tilt 
channel.  The  exponential  hold  in  Fig.  4.6-1  has  continous 
dynamics  identical  to  those  of  the  tilt  disturbance;  its 
purpose  is  to  smooth  the  sampled  z between  updates.  The 
tilt  mirror  and  feedback  loop  are  identical  for  the  two 
systems;  the  control  specifications  and  control  design  are 
described  in  Section  4.4.  The  scalar  gain  K provides  a unity 
DC  gain  from  C^t^)  to  in  the  classical  system. 


Figure  4.6-1  Optimal  Model  Following  Controller 

for  Turbulence-Induced  Tilt 


DISTURBANCE  l 

II nnc.  CONTROl  SYSTEM 
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Figure  4.6-3  shows  the  nulling  error  variances  for 
the  two  systems  as  a function  of  control  sample  rate.  The 
autocorrelation  function  for  a first  order  shaping  filter 
(9  Hz  bandwidth)  was  used  to  derive  the  third  curve;  this 
curve  establishes  an  upper  bound  on  transport  lag  effects. 

The  continuous  system  lines  are  asymptotic  limits  of  the 
digital  system  curves  for  high  sample  rates. 

The  optimal  system  shows  consistently  superior  per- 
formance at  all  sample  rates  and  exceeds  the  performance  of 
the  continuous  classical  system  at  a sample  frequency  cf  1200  Hz 
The  optimal  system  offers  significant  advantages  over  the  classi 
cal  in  data  processing  (sample  rate)  requirements. 

In  addition  to  an  alleviation  of  sample  rate  require- 
ments the  optimal  controller  can  maintain  good  performance  with 
lower  grade  hardware  (lower  response  bandwidth).  Figure  4.6-4 


NULLING  ERROR  VARIANCE  - TURBULENCE  VARIANCE  RATIO  FOR  TILT  CHANNEL 

wrius 

CONTROLLER  SAMPLE  RATE 


I 


Figure  4.6-3 


Nulling  Error  Comparison  Between  Classical 
and  Optimal  Controllers 


TRACKING  ERROR  VARIANCE  - TURBULENCE  VARIANCE  RATIO  FOR  TILT  CHANNEL 


Figure  4.6-4  Nulling  Errcr  Variance  Versus  Tilt 

Mirror  Bandwidth 

shows  the  nulling  error  variance  of  the  tilt  channel  controller 
as  a function  of  tilt  mirror  bandwidth.  The  error  variance 
reaches  a point  of  diminishing  return  at  500  Hz;  there  is  little 
to  be  gained  by  increasing  mirror  bandwidth.  Since  a higher  re- 
sponse bandwidth  generally  implies  high  quality  (expensive)  hard- 
ware, use  of  an  optimal  controller  could  have  a significant  impact 
on  component  development  goals  and  cost. 

The  optimal  system  is  quite  insensitive  to  mismatches 
between  the  design  disturbance  model  (used  to  design  controller 
and  holding  circuit)  and  the  truth  model.  Figure  4.6-5  shows 
that  design  model  bandwidth  errors  of  z20ao  yield  insignificant 
changes  in  nulling  error  variance.  The  bandwidths  of  the  tur- 
bulence disturbances  are  proportional  to  the  local  wind  velocity 
(Ref.  1);  hence,  the  optimal  controller  might  be  thought  to  be 
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TRACKING  ERROR  - TUflflULENCE  STANDARD  DEVIATION  RATIO  FOR  TILT  CHANNEL 
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Figure  4.6-5  Tilt  Channel  Mismatch  Sensitivity 

tolerant  to  r.20%  errors  in  measured  wind  velocity;  however, 
errors  in  the  turbulence  model  bandwidths  will  also  be  present 
due  to  uncertainty  in  the  modeled  relationship  between  wir..d 
velocity  and  turbulence  bandwidth. 

In  summary,  the  tilt  channel  optimal  controller  offers 


• Disturbance  nulling  performance  superior 
to  that  of  a typical  classical  design 

• Compensation  of  low  bandwidth  hardware 
to  yield  performance  comparable  to  that 
of  higher  quality  components 

• Insensitivity  of  the  design  to  disturbance 
model  mismatches. 


4.6.2  Deformable  Mirror  Control  Evaluation 


In  early  evaluation  of  the  deformable  mirror  it  was 
found  that  a final  scalar  adjustment  of  the  gain  matrix  derived 
from  the  pseudo-optimal  control  results  (summarized  in  Section  4. 
is  required  to  achieve  proper  disturbance  nulling.  The  inade- 
quacy of  the  pseudo-optimal  feedforward  control  is  evidenced 
in  Fig.  4.6-6  which  shows  an  astigmatism  time  history  and 
the  mirror  output  (negated)  that  attempts  to  counteract  it. 

The  mirror  output  is,  very  nearly,  a scalar  attenuation  of 
the  desired  signal;  this  attenuation  is  a result  of  a spurious 
restriction  on  control  magnitude  imposed  by  the  R”1  matrix  of 
the  fictitious  regulator.  Consider  the  relationship  for  the 
diagonal  elements  of  the  R-1  matrix,  i.e.,  those  elements  which 
determine  the  magnitude  of  the  primary  disturbance  nulling  con- 
trol commands. 


r|  55  RiJ  = “ ■ C2(  1-S2  ) (4.5-1) 

The  magnitude  of  this  terra  determines  how  much  of  the  distur- 
bance states  "gets  through"  to  command  the  mirror.  For  example, 
with : 


C - .1 

r*  =.  --Q396 
1 qi 

The  maximum  value  of  r*  occurs  when 

for  ; «*  /2/2 


(4.6-1) 

(4.6-2) 


✓'2  / 2 ; i.e.: 


(4.6-3) 


If  the  damping  ratio  of  the  mirror  is  less  than  Z2/2, 


Figure  4.6-6  Disturbance  and  Deformable  Mirror  Time  Histories 


the  fictitious  regulator  derived  will  represent  an  optimization 
problem  for  which  there  is  a restriction  on  control  magnitude; 
i.e.,  the  value  of  r“  serves  to  attentuate  control  commands. 


A simple  artifice  to  relax  this  undesirable  restriction 
is  to  scale  up  the  derived  R-1  matrix  by  the  factor: 


S.F. 


-i 


(4.6-4) 


This  scaling  has  the  effect  of  making  the  diagonal  elements 
of  the  R-^  matrix  equal  to  those  for  the  i * v'2’/2  case  (which 
imposes  no  attenuation  of  control  commands)  and  increasing  the 
off-diagonal  elements  to  maintain  the  proper  proportion  of 
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crossfeed.  The  success  of  the  R-^  scaling  artifice  is  illus- 
trated by  the  results  of  the  following  example  problem. 

Figures  4.6-7  and  4.6-8  are  block  diagrams  which  show, 
respectively,  the  optimal  control  system  for  the  deformable  mir- 
ror developed  in  Section  4.5  and  a classical  inverse-filter 
control  system. 


Figure  4.6-7  Deformable  Mirror  Optimal  Control  System 


TURBULENCE  MOOElS 


I 

INVERSE  FILTER  CONTROLLER  | MIRROR  OYNAMiCS 

I 

I 


Figure  4.6-8  Classical  Inverse  Filter  Control  System 
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The  inverse  filter  computes  controls  assuming  that  the  defor- 
mable mirror  reaches  steady  state  between  samples:  i.e.  with 
the  mirror  dynamics  defined  by : 

x = Fx  + Gu  (4.6-5) 

in  steady  state: 

0 « Fx(tk)  + Gu(tk)  (4.6-6) 

and 

*<V'  - (4-6-7) 


therefore : 


u(tk)  ■-  (GTG)-1GTF(HinHmT)-1Hm£(t,t) 

(4.S-8) 

The  inverse  filter  approach  depends  primarily  on  high 
mirror  bandwidth  to  achieve  the  steady  state  assumption  and 
high  sample  rate  to  avoid  transport  lag  effects. 

Primary  control  of  the  two  astigmatism  disturbances 
with  suppression  of  mirror  focus  cross-coupling  was  used  as 
a test  case  to  compare  the  performance  of  the  two  systems. 

The  astigmatism  truth  models  chosen  were  those  for  the  phase 
conjugate  scenario.  For  very  low  distrubance  bandwidths , 
such  as  those  of  the  AMD  scenario,  the  optimal  controller 
for  the  3300  Hz  bandwidth  deformable  mirror  reduces  (very 
nearly)  to  an  inverse  filter  controller;  the  classical  and 
optimal  controllers  are  identical.  Amplitude  cross-couplings 
into  focus  of  2°o  and  3^,  respectively,  for  the  two  astigmatism 
modes  were  included  in  the  mirror  dynamics;  i.e.  the  coefficient 
coupling  amplitude  of  optical  mode  4 into  optical  mode  3 is 


I 

o 2 

(-.02)oj  ; the  coefficient  for  mode  5 into  mode  3 is  (-.03)u>m. 
m m 

Figures  4.6-9  through  4.6-11  summarize  the  error 
nulling  performance  of  the  two  systems.  Figure  4.6-9  shows 
the  nulling  error  variances  of  the  two  astigmatism  modes  as 
a function  of  damping  ratio;  the  optimal  controller  out -per- 
forms the  inverse  filter  controller  at  all  damping  ratios. 
Figure  4.6-10  shows  the  variance  of  the  focus  mode  produced 
by  cross-coupling.  The  curve  for  the  optimal  controller  shows 

a very  low  variance  for  the  focus  coupling:  roughly  10~  m 

-12  ^ 

as  opposed  to  10  ra"1  for  atmospheric  turbulence  focus  distur- 
bances. The  focus  cross-coupling  variances  for  the  inverse 
filter  controller  are  roughly  two  orders  of  magnitude  lower 
than  those  of  the  optimal  system;  the  primary  control  responses 
of  the  optimal  controller  produce  cross-couplings  which  are 
more  difficult  to  suppress  than  those  produced  by  the  piece- 
wise  constant  inverse  filter  commands.  Figure  4.6-11  shows 
the  normalized  intensity  degradation  for  the  two  systems, 
i.e.  the  intensity  degradation  produced  with  control  of  as- 
tigmatism divided  by  that  produced  with  no  control  of  astig- 
matism. The  optimal  controller  again  outperforms  the  inverse 
filter;  .052  normalized  intensity  degradation  for  the  optimal 
system  versus  .072  for  the  inverse  filter. 

In  summary,  the  optimal  controller  for  the  deformable 
mirror  has  the  following  qualities: 

• Primary  control  of  astigmatism  superior  to  that 
of  an  inverse  filter  controller 

• Adequate  suppression  of  cross-coupling  dynamics 

i 

j 

• Lower  overall  intensity  degradation  chan  inverse 
filter  controller. 

t 

• i 
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MIRROR  DAMPING  RATIO 


Figure  4.6-11  Normalized  Intensity  Degradation 


5. 


INTEGRATED  ESTIMATION-CONTROL  SIMULATION  RESULTS 


In  this  chapter  results  of  the  performance  evaluation 
of  the  integrated  estimation-control  algorithms  for  AMD  and  PC 
are  presented.  In  addition  to  full  disturbance  estimation 
(five  truth  model  optical  modes  and  five  estimator  modes),  a 
number  of  limited  mode  estimation  and  control  cases  were  run. 
For  example,  the  PC  system  was  implemented  with  estimation  and 
control  of  the  two  tilt  modes  and  focus  while  the  truth  model 
included  all  five  optical  disturbance  modes.  An  interim  iiver- 
gence  of  the  AMD  estimator  was  discovered  during  the  integrated 
simulations;  the  source  of  this  divergence  and  methods  to  avoid 
it  are  discussed  in  this  chapter,  also. 


5.1  INTEGRATED  PHASE  CONJUGATE  SYGTEM  RESULTS 

Figure  5.1-1  shows  the  integrated  estimation  and 
control  simulation  results  for  the  phase  conjugate  system  as 
a function  of  sensor  noise  level.  In  addition  to  the  curve 
shown  for  the  five  mode  case,  curves  are  provided  for  simula- 
tions in  which  only  modes  1,  2,  and  3 (xtilt,  ytilt,  focus) 
and  only  modes  1 and  2 were  controlled.  In  all  cases  the 
truth  model  included  all  five  aberration  modes.  A curve  for 
the  Kalman  filter  with  perfect  control  and  the  no  estimation 
and  control  line  are  included  for  comparison. 

The  five  mode  estimation-control  curve  shows  a 5 to 
7 percent  intensity  degradation  from  the  Kalman,  filter  with 
perfect  control  line;  this  is  a reasonable  degree  of  degrada- 
tion considering  the  results  of  Section  4.6.1.  That  is, 
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Figure  5.1-1  Phase  Conjugate  System  Performance 


roughly  1 to  2%  nulling  error  degradation  is  to  be  expected 
for  each  of  the  tilt  channels;  similar  degradations  should  be 
expected  for  the  focus  and  deformable  mirror  channels.  As  a 
consequece  of  the  separation  theorem  (Ref.  2),  the  five  mode 
estimation  and  control  results  represe.it  the  performance  of 
the  optimal  system  consistent  with  the  modeling  assumptions 
and  parameters  used  in  the  present  study. 

Additional  simulations  were  run  to  evaluate  the 
system  performance  for  limited  disturbance  estimation  and 
control.  Again,  all  five  disturbance  models  were  present  in 
the  truth  model;  the  filtering  algorithm  was  structured  to 


estimate  only  limited  combinations  of  aberration  modes.  No 
alteration  to  the  filtering  algorithm  was  made  to  account  for 
the  presence  of  the  neglected  disturbances;  it  was  assumed 
that  the  orthogonality  of  the  aberration  modes  and  the  spatial 
filtering  properties  of  the  shearing  interferometer  itself 
would  "wash  out"  undesired  errors  due  to  neglected  higher 
order  aberrations.  The  mode  1,  2,  3 case  in  which  astigma- 
tism was  neglected  shows  Strehl  ratio  degradations  from  the 
five  mode  correction  case  ranging  from  1 percent  at  the  high 
sensor  noise  value  to  4 percent  at  the  low  sensor  noise  level. 
The  case  in  which  only  the  two  tilt  modes  were  controlled 
shows  degradations  of  2 percent  at  high  noise  to  6 percent  at 
low  noise  levels.  One  can  infer  from  these  results  that  at 
high  sensor  noise  levels  there  i$  little  to  be  gained  by  esti- 
mation and  control  of  higher  order  aberration  modes;  indeed, 
the  savings  in  required  hardware  (i.e.,  focus  and  deformable 
mirrors)  may  justify  neglect  of  these  modes. 


5.2 


INTEGRATED  ACTIVE  MULTIDITHER  SYSTEM 


5.2.1  Test  Cases  and  Simulation  Results 

Integrated  estimation  and  control  simulations  were 
performed  for  the  AMD  system  of  a dither  frequency  of  250  Hz. 

The  500  Hz  dither  was  not  evaluated  for  two  reasons: 

• The  tilt  and  focus  mirror  bandwidths  specified 
(1000  Hz  and  500  Hz,  respectively)  would  yield 
poor  following  of  a 500  Hz  dither  signal. 

• Preliminary  estimation  evaluation  (Section  3) 
indicated  that  with  speckle  performance  was 
nearly  identical  at  250  Hz  and  500  Hz  dither 
frequencies. 

The  first  concern  could  be  obviated  by  dither  command  signal 
shaping  (i.e.  a "lead'1  prefilter),  increased  dither  amplitude, 
or  by  implementation  of  a separate  high  bandwidth  dither  mirror. 
The  dither  signal  processing  approaches  would  require  experimen- 
tation; the  hardware  alternative  would  require  a sophisticated 
and  expensive  optical  component.  The  estimation  performance 
with  speckle,  the  most  realistic  representation  of  the  AMD 
system,  indicates  negiigable  improvement  at  the  higher  dither 
frequency,  efforts  to  dither  at  500  Hz  will  yield  little,  if 
any,  reward. 

The  performance  of  the  integrated  AMD  system  at  a 
range  of  500  km  is  indicated  in  Fig.  5.2-1.  The  no  speckle 
case  and  the  no  estimation  and  control  case  establish  upper 
and  lower  bounds,  respectively,  on  system  performance.  In 
the  presence  of  speckle  the  integrated  estimation  and  control 
system  achieves  a Strehl  ratio  of  .550:  a significant  improve- 
ment over  the  no  estimation  and  control  value  of  .300.  but  not 
a very  encouraging  result  from  the  standpoint  of  overall  system 
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Figure  5.2-1  Integrated  AMD  Performance, 

Range  * 500  km 


efficiency.  Again,  speckle  effects  cause  a very  severe 
degradation  in  measurement  -quality , reducing  the  effective 
signal-to-noise  ratio  of  the  photodetector  to  roughly  3.  With 
measurement  noise  'levels  of  such  high  magnitude  even  an  optimal 
system  cannot  offer  totally  favorable  results.  It  should  be 
mentioned  that  the  estimator  underwent  an  interim  divergence 
in  the  no  speckle  case  (no  divergence  was  detected  with 
speckle).  This  divergence  was  avoided  by  reinitialization  of 
the  filtering  algorithm;  the  Strehl  ratio  shown  in  Fig.  5.2-1 
corresponds  to  the  "fixed"  filter.  The  nature  and  correction  of 
the  filter  divergence  are  discussed  in  Section  5.2.2. 

Figure  5.2.2  shows  the  performance  of  the  AMD  system 
at  100  km  range.  Again,  the  no  speckle  and  no  control  results 
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Figure  5.2-2 


Integrated  AMD  Performance, 
Range  * 100  km 


are  presented  for  comparison;  the  no  speckle  results  are  con- 
siderably better  than  those  for  the  500  km  range  as  a result 
of  the  smaller  transit  delay  lag  and  an  increased  signal-to- 
noise  ratio  at  the  closer  range.  As  the  estimation  and  con- 
trol with  speckle  line  indicates,  speckle  washes  out  any  im- 
provement in  the  signal-to-noise  ratio,  thereby  degrading 
performance  to  a Strehl  ratio  of  .63.  No  divergence  of  the 
filter  was  detected  in  either  the  no  speckle  or  speckle  cases 

5.2.2  Filter  Divergence 


Filter  divergence  (instability)  is  a common  problem 
of  extended  Kalman  filters  and  generally  arises  from  the 
linear  approximations  used  in  the  filtering  algorithm.  In 
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the  AMD  estimator  there  is  a linearization  of  the  measurement; 
this  approximation  is  an  additional  source  of  error  to  the 
filter  estimates  that  is  not  accounted  for  explicitly  in  the 
algorithm.  A divergence  of  the  estimation  algorithm  occurred 
in  the  500  km  range  with  no  speckle  simulation.  While  the 
no  speckle  simulations  do  not  represent  physically  realistic 
conditions,  the  nature  and  correction  of  the  divergence  should 
be  examined  to  determine  why  no  divergence  occurred  with  speckle 


The  time  history  of  the  turbulence  mode  strengths 
and  the  corresponding  filter  estimates  near  the  divergence  is 
shown  in  Fig.  5.2-3.  The  mode  strengths  and  estimates  all  ap- 
proach a small  value  simultaneously;  as  the  mode  strengths  move 
apart  the  estimates  begin  to  oscillate  with  increasing  amplitude 
(diverge).  A similar  problem  was  observed  by  Lewantowicz  (Ref.  12) 
in  which  a conical  scan  estimator  developed  a sign  error  in  the 
estimate  and  diverged. 


Figure  5.2-3(a)  Divergent  Filter  Time  History 
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There  are  a number  of  artifices  that  can  be  used  to 
correct  divergence  of  an  extended  Kalman  filter.  In  Ref.  11, 
for  example,  the  filter  residual  was  monitored;  if  the  resi- 
dual appeared  to  be  growing,  a sign  error  would  be  assumed 
and  the  sign  of  the  estimate  changed.  A similar  logic  could 
be  applied  to  the  present  case;  the  residual  could  be  monitored 
to  detect  the  divergence  and  the  filter  reinitialized.  Table 
5.2-1  shows  the  estimation  errors  for  the  divergent  and  "fixed" 
filters  (dynamic  logic  was  not  implemented;  the  filter  was 
simply  reinitialized  at  the  time  the  divergence  occurred). 

The  estimation  error  for  the  focus  mode  (Mode  3)  is  still 
quite  large;  this  is  a result  of  a very  large  defocus  excur- 
sion (approximately  10  peak  with  a time  constant  of  0.15 
sec)  following  reinitialization.  The  estimation  errors  shown 
in  Table  5.2-1  imply  that  no  additional  divergence  of  the  fil- 
ter occurred  after  reinitialization. 


TABLE  5.2-1 

ESTIMATION  ERRORS  FOR  DIVERGENT 
AND  FIXED  FILTERS 


MODE 

1 

2 

3 

4 

5 

Disturbance 

rms  (cd)  (um) 

1.20 

1.20 

0.23 

0.33 

0.33 

Estimation 

Divergent 

1.07 

1.00 

0.22 

0.30 

0.30 

Errors  (um) 

Reinitialized 

0.69 

0.46 

0.30 

0.17 

0.22 
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No  divergence  was  detected  in  simulations  which  in- 
cluded speckle.  As  outlined  in  Section  3.1,  the  filter  was 
altered  to  account  for  speckle  by  making  a sizeable  increase 
in  the  assumed  detector  noise;  this  increase  causes  the  filter- 
ing algorithm  to  put  less  weight  on  the  measurements  relative 
to  the  extrapolations  in  computing  the  estimates.  This  approach 
to  dealing  with  speckle  introduces  additional  stability  into 
the  filtering  algorithm;  this  increased  stability  also  ac- 
commodates the  errors  resulting  from  the  linearization  of  the 
measurement . 


6. 


BEAM  CLEANUP  SYSTEM 


The  beam  cleanup  controller,  used  to  control  phase- 
front  aberrations  generated  by  the  laser  device  itselt,  is 
described  in  this  chapter.  The  design  discussed  here  deals 
only  with  the  control  problem;  since  the  laser  power  levels 
involved  exceed  the  sensor  noise  levels  by  several  orders  of 
magnitude,  "perfect"  state  reconstruction  is  possible  via  a 
Luenberger  state  estimator  (observer,  Red.  2).  The  .first 
section  of  this  chapter  will  deal  with  the  development  of  an 
optimum  error  budget;  in  the  second  section  a feedforward 
controller-deformable  mirror  system  is  designed  to  meet  the 
optimum  error  budget  in  a best -least -squares  sense;  simulation 
results  and  evaluation  of  the  control  system  are  presented  in 
the  third  section. 


6.1  OPTIMUM  ERROR  BUDGET 
• 

The  prime  limitation  of  a beam  cleanup  control  system 
stems  from  the  fact  that  a single  deformable  mirror  must  at- 
tempt to  null  phasefront  aberrations  associated  with  a number 
of  different  radiation  wavelengths.  In  the  present  case,  for 
example,  tho  total  laser  power  is  divided  among  three  distinct 
wavelengths  or  radiation  spectrum  lines;  the  single  deformable 
mirror  cannot  null  simultaneously  the  aberrations  of  all  three 
lines . 


Since  perfect 
best  alternative  is  to 
that  will  minimize  the 


disturbance  nulling  is 
find  a rime  history  of 
intensity  degradations 


impossible,  the 
mirror  output 
due  to  the 
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disturbances,  Eq.  6.1-1  is  a quadratic  performance  index  which 
is  appropriate  to  this  design  goal: 


j U£1+c  ) Q1(£1+£ )+(£2+c  ) Q2(c2+c  )+(c3+c)  Q3(c_3+c  )]  dt 
° (6.1-1) 


or,  in  more  compact  notation: 


s 


3 

Z 

i*l 


(^i 


^)TQi(-l 


+ c ) dt  (6.1-2) 


Here  c,  is  the  vector  of  optical  mode  disturbances 
associated  with  the  i radiation  spectrum  line,  c is  the 
vector  of  deformable  mirror  outputs  which  attempts  to  counter- 
act the  phasefront  aberrations  of  the  three  lines.  The  s>i's 
are  the  intensity  degradation  matrices  derived  in  Section 
4.3.1;  these  are  diagonal  matrices  which  are  functions  of 
both  radiation  power  and  wavelength  of  a given  line;  i.e. 
for  the  diagonal  elements  are: 


Q » J Z. 

qll  max,  4 

q22  * ^max^  4 
k2 

q33  * Imaxi 


(4.3-9) 
(4.3-10) 
(4 . 3-11) 


k 


i 


(2.3-2) 


The  problem  is  to  find  the 
minimizes  the  performance  index  Eq . 
J with  respect  to  5 one  obtains: 


mirror  output,  6,  that 
6.1-1.  Differentiating 


0 J 
3c 


dt 


(6.1-3) 
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The  min  imam  degradation  occurs  when  ^ is  equal  to  a zero 
vector;  this  implies  that  the  integrand  of  Eq . 6.1-3  is  zero 
for  all  t : 


» 


i 


i 

< i 

'! 

V 


3 

l Q^c,  + Qic  - 0 (6.1-4) 

i"l 


3 3 

C l Qt)c  --  l Qtc 

i»l  i=»l  1 

3 

3 

(6.1-5) 

/.  c * - ( l Q.,) 
i-1 

I,  Qi-i 

i”l 

(6.1-6) 

Eq.  6.1-6  defines  the  optimum  time  history  of  the 
deformable  mirror  output  in  counteracting  the  phasefront  aber- 
rations of  the  three  lines.  Since  the  Q^'s  are  all  equal  up  to 
a multiplicative  factor,  further  simplification  of  Eq . 6.1-6  is 
possible , i . e . ; 


Q±  * 

I 

max 

.kiQ 

(6.1-7) 
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(6.1-8) 
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3 

l 

i*l 

Qi 

- PtQ 
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Pt  * 

3 

l 
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(6.1-12) 
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Substituting  Eq . 6.1-9  and  6.1-10  into  Eq . 6.1-6  one  obtains: 


£ - - l ~ c . (6.1-13) 

i-1  t 

Eq . 6.1-13  indicates  that  the  deformable  mirror  should  follow 
the  power  weighted  sum  of  the  phasefront  aberrations  of  the 
three  lines. 


6.2  THE  BEAU  CLEANUP  CONTROL  SYSTEM/ SIMULATION  MODEL 

A control  system  will  now  be  designed  to  follow  the 
optimum  output  time  history  derived  in  Section  6.1.  The  de- 
sign will  be  a direct  adaptation  of  the  feedforward  controller- 
deformable  mirror  system  derived  in  Section  4.5.  The  design 
approach  is  to  derive  the  optimal  feedforward  controller 
for  the  phasefront  disturbances  of  each  line,  by  itself;  the 
single  line  controller  outputs  are  then  attenuated  by  the  frac- 
tion of  the  total  laser  power  of  that  line.  Figure  6.2-1  is  a 
block  diagram  of  the  resulting  control  system. 

The  control  system  as  shown  is  analog.  The  disturbance 
bandwidths  involved  were  in  the  kilohertz  range;  a digital  con- 
trol implementation  would  require  prohibitively  high  sample  rates. 
Again,  the  phasefront  aberration  measurements  are  assumed  to 
be  perfect;  in  steady  state  the  x^'s  used  in  computing  the 
control  will  be  identical  to  the  states  of  the  disturbance 
model . In  the  computer  simulation  of  the  beam  cleanup  system 
the  truth  model  disturbance  states  are  used  directly  in  the 
control  computation;  thus,  the  simulation  results  of  Section 

6.3.2  will  reflect  the  steady  state  performance  cf  the  system. 
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6.3 


BEAM  CLEANUP  SIMULATION 


6.3.1  Test  Case  Parameters 


The  test  case  used  to  evaluate  the  beam  cleanup 
design  is  described  in  this  section.  The  problem  is  to  con- 
trol the  phasefront  aberrations  of  three  radiation  spectrum 
lines;  the  wavelengths  and  radiation  power  levels  are  given 
in  Table  6.3-1. 


« 


TABLE  6.3-1 

RADIATION  WAVELENGTHS  AND  POWER  LEVELS 

T-1848 


LINE  NO. 

WAVELENGTH  (urn) 

P.  /P 
i'  t 

1 

3.714 

. 113 

2 

3 . 800 

. 792 

3 

3.927 

. 085 

Phasefront  disturbances  for  each  line  are  modeled 
by  the  dynamics  of  three  optical  modes  (two  tilt  modes  and 
focus).  The  spectrum  for  a typical  phasefront  disturbance 
mode  is  shown  in  Fig.  6.3-1;  break  frequencies  and  rms  levels 
for  all  nine  disturbance  models  are  listed  in  Table  6.3-2. 
These  spectra  are  simulated  in  the  time  domain  by  the  outputs 
of  linear  shaping  filters  driven  by  white  noise. 

6.3.2  Simulation  Results 


Figure  6.3-2  shows  the  efficiency  of  the  beam  cleanup 
system  as  a function  of  deformable  mirror  bandwidth.  The 
asymptote  labeled  "perfect  control"  indicates  that  the  rms 
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Figure  6.3-1  Phasefront  Disturbance  Spectrum 


TABLE  6.3-2 

DISTURBANCE  SPECTRA  PARAMETERS 


LINE  NO. 


MODE 


T-1849 


RilS  (um) 


XTILT 

2tt(1000) 

YTILT 

2tt(  1000 ) 

FOCUS 

2 it  ( 1200  ) | 

XTILT 

2^(800) 

YTILT  | 2t( 800 ) 


FOCUS  | 2t( 900 ) 


XTILT  | 2tt(1200) 


YTILT  2* (1200) 


FOCUS  2*t  ( 1750  ) 


2tt  ( 1400  ) 2-t(15) 


2tt  ( 1400 ) 2t(15) 


2^(1500)  j 2^(20) 


2tt  ( 1000  ) ! 2^(60) 


2^(1000)  ! 2t(60) 


2*r  ( 1100  ) | 2-r(40) 


2t  ( 1450 ) | 2*t  ( 12  ) 


2^ ( 1450 ' 


2t  ( 12) 


2 tt  ( 2000  ) j 2- (50) 


1 
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Figure  6.3-2  Beam  Cleanup  System  Performance 

levels  chosen  cause  a severe  intensity  degradation;  even 
with  exact  following  of  the  optimal  time  history  derived  in 
•Section  6.1  roughly  55%  of  the  potential  target  intensity  is 
lost.  Although  the  disturbance  levels  may  be  unrealistically 
high,  the  example  problem  does  give  some  indication  of  hardware 
requirements  for  the  beam  cleanup  system.  With  a mirror  band- 
width of  3300  Hz,  the  bandwidth  assumed  for  the  PC  and  AMD 
systems,  the  target  intensity  is  degraded  to  .225;  even  with 
an  optimal  controller  the  mirror  is  simply  too  slow  to  follow 
| the  high  bandwidth  disturbances.  Indeed,  a mirror  bandwidth 

near  10  kHz  would  be  required  for  disturbance  nulling  approach- 
i fn8  the  maximum  possible  at  the  specified  disturbance  bandwidth. 

' t 
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The  example  problem  results  also  indicate  that  the 
performance  index,  Eq.  6.1-1,  is  inappropriate  in  situations 
having  high  rms  levels  of  phasefront  distortion.  Equation 
6.1-1  is  based  on  a quadratic  approximation  of  the  intensity 
pattern  at  the  target  (see  Section  4.3.1);  this  approximation 
is  accurate  only  near  the  origin  of  the  intensity  function. 

For  high  levels  of  phasefront  aberration  the  quadratic  approxi- 
mation breaks  down  and  may  yield  normalized  intensity  degrada- 
tions which  exceed  unity  as  indicated  by  Fig.  6.3-3. 


Figure  6.3-3  Gaussian  Intensity  Pattern  and 

Quadratic  Approximation 


The  inaccuracy  of  the  quadratic  approximation  ad- 
versely affects  the  performance  index  and  the  optimal  control 
time  history  derived  from  it.  Two  approaches  to  the  deriva- 
tion of  a more  accurate  control  methodology  are: 


Use  the  assumed  Gaussian  intensity 
function  (Eq.  2.3-1)  to  characterize 
target  intensity  degradations  in  the 
performance  index,  i.e.: 


• Null  only  the  phasefront  disturbances  of 
the  line  having  the  largest  fraction  of 
the  total  power. 

The  first  alternative  is  the  more  precise  of  the  two; 
the  [~Imax  -g.(c. +5)1  terms  represent  a more  accurate  repre- 
sentation  of  intensity  degradations  than  the  quadratic  approxi- 
mations used  in  Eq . 6.1-1.  Since  the  performance  index  is  no 
longer  quadratic,  the  minimization  of  J is  no  longer  a simple 
matter  of  taking  the  derivative  of  the  performance  index  with 
respect  to  £ and  setting  it  equal  to  zero.  A time  invariant 
error  nulling  budget  similar  to  Eq.  6.1-13  could  be  derived 
by  parameter  optimization  methods;  i.e.,  the  optimal  control 
would  be  a weighted  sum  of  the  phasefront  aberrations  of  the 
three  lines.  The  weighting  factors  would  have  to  be  determined 
through  monte-carlo  analysis;  this  can  be  a very  costly  and 
time  consuming  process.  While  admittedly  suboptimal,  the 
second  approach  is  simple  to  derive,  implement,  ar.d  analyze. 

In  situations  similar  to  the  test  case,  in  which  80^  of  the 
laser  energy  was  in  a single  spectral  line,  adequate  perfor- 
mance can  be  obtained.  With  a mirror  and  control  system  capa- 
ble of  accurately  nulling  the  phasefront  aberrations  of  the 
high  power  spectral  line,  nearly  80*^  of  the  laser  power  could 
be  rendered  unaberrated  (i.e.,  nearly  all  of  the  power  in  the 
3.8  urn  line  could  be  cleaned  up).  A normalized  maximum  in- 
tensity of  0.8  is  a considerable  improvement  over  the  value 
of  0.45  obtained  with  the  controller  based  on  Eq.  6.1-13. 
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6.4  CHAPTER  SUMMARY 


The  beam  cleanup  design  and  evaluation  indicates: 


• For  low  disturbance  levels  (i.e.,  levels 
for  which  the  quadratic  approximation 
(Eq.  6.1-1)  is  accurate)  the  optimal  beam 
cleanup  time  history  is  the  power  weighted 
sum  of  the  phasefront  disturbances  of  the 
three  lines. 

• A mirror  bandwidth  of  10  kHz  is  required 
for  accurate  optimal  time  history  tracking 
of  disturbances  in  the  1 kHz  range. 

• With  high  disturbance  levels  (i.e.,  levels 
for  which  the  quadratic  approximation 
(Eq.  6.1-1)  is  not  accurate)  the  optimal 
beam  cleanup  time  history  is  no  longer  the 
power  weighted  sum  of  the  phasefront  dis- 
turbances of  the  three  lines.  The  optimal 
weighting  of  the  phasefront  disturbances  of 
the  three  lines  must  be  derived  through 
monte  carlo  analysis  using  a Gaussian  per- 
formance index. 


i 


l 

1 


7. 


CONCLUSIONS  AND  RECOMMENDATIONS 


i 


7.1  CONCLUSIONS 

The  extended  Kalman  filter  developed  in  the  present 
program  is  an  effective  algorithm  for  active  multidither  phase- 
front  estimation.  The  Kalman  filter  successfully  estimates 
the  states  of  five  turbulence  induced  phasefront  disturbance 
modes  in  the  presence  of  severe  speckle  and  beam  transit  de- 
lays. The  performance  of  the  extended  Kalman  filter  (in  terms 
of  an  average  Strehl  ratio  for  the  scenario)  was  found  to  be 
essentially  the  same  for  dither  frequencies  of  250  Hz  and  500 
Hz;  the  use  of  the  lower  dither  frequency  reduces  the  mirror 
bandwidth  required  to  follow  the  dither  commands, 

The  Kalmar  filter  developed  for  phase  conjugate  aber- 
ration estimation  yielded  Strehl  ratio  increases  of  0.05  to 
0.28  (depending  on  sensor  noise)  over  a candidate  least-squares- 
fit  low  pass  filter  estimator.  Sensitivity  studies  to  deter- 
mine the  effects  of  truth  model/filter  signal  to  noise  mis- 
matches which  may  result  from  target  heating  indicate  negli- 
gible degradation  of  estimator  performance. 

Optimal  disturbance  nulling  control  was  shown  to 
yield  significant  performance  improvements  over  typical  clas- 
sical designs.  Superior  performance  of  the  optimal  controller 
is  attributable  to  the  use  of  the  full  disturbance  state  (i.e., 
states  related  to  the  time  derivatives  of  the  disturbance  pro- 
cess as  well  as  the  disturbance  magnitude)  in  computing  control 
commands.  Using  the  time  derivative  states,  the  controller 
can  compute  commands  which  anticipate  the  future  of  the  dis- 
turbance process  in  addition  to  nulling  present  disturbances. 

The  optimal  controller  for  the  tilt  channel  outperformed  a 
candidate  classical  controller  at  all  sample  rates  and  had 
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nulling  error  variances  lower  than  that  of  a continuous  class- 
cal  controller  for  sample  rates  higher  than  1.2  kHz.  Similarly, 
an  optimal  feedforward  controller  designed  for  the  deformable 

1 mirror  outperformed  a classical  inverse-filter  controller  in 

nulling  astigmatism  errors  while  suppresing  mirror  surface 
figure  cross-coupling.  Parametric  studies  showed  that  use  of 
optimal  control  algorithms  can  relax  hardware  performance  re- 
quirements; for  example,  increases  in  tilt  mirror  bandwidth  j 

beyond  500  Hz  yield  negligible  improvement  is  distrubance  nulling.  j 

Optimal  controller  performance  also  was  found  to  be  insensi- 
tive to  20  percent  mismatches  of  the  design  model  and  truth 
model  disturbance  bandwidths. 

Simulations  of  the  integrated  estimation  and  control 
systems  indicate  good  performance:  non-ideal  control  had  negli- 
gible effect  on  convergence  of  the  active  multidither  estimator; 
non-ideal  control  reduces  the  Strehl  ratio  produced  by  the  phase 
conjugate  system  0.05  to  0.07  (depending  on  sensor  noise)  be- 
low the  corresponding  ideal  control  results. 

The  results  of  the  beam  cleanup  system  design  and 
evaluation  indicated  that  the  optimal  control  time  history  de- 
pends on  the  phasefront  disturbance  levels  as  well  as  the 
power  contained  within  each  line.  For  low  disturbance  levels,  t 

i.e.,  levels  for  which  the  Gaussian  intensity  can  be  approxi- 
mated by  a quadratic  function,  the  optimal  beam  cleanup  time 
history  is  the  power  weighted  sum  of  the  phasefront  distur- 
bances of  the  three  spectral  lines.  At  disturbance  levels  t 

for  which  a quadratic  approximation  is  inaccurate,  the  optimal 
weighting  of  the  phasefront  disturbances  must  be  derived  * I 

J through  monte  carlo  analysis  using  a Gaussian  performance  in-  f 

dex.  Parametric  studies  also  indicated  that  a mirror  band- 
I width  of  10  kHz  is  required  for  accurate  optimal  time  history  < 

i f tracking  of  disturbances  in  the  1 kHz  range. 
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In  summary,  optimal  estimation  and  control  have  been 
shown  to  offer: 

• An  effective  approach  to  active  multidither 
estimation  in  the  presence  of  severe  speckle 

• Phase  conjugate  estimation  superior  to  classi- 
cal methods 

• Disturbance  nulling  control  performance  superior 
to  that  of  typical  classical  designs 

• Insensitivity  of  control  performance  to  distur- 
bance model  mismatches 

• Potential  relaxation  of  active  component  band- 
width requirements 

• Spatial,  temporal  optimization  of  deformable 
mirror  control 


7 . 2 RECOMMENDAT IONS 

Isoplanaticism  Correction  - At  long  ranges  (i.e., 

100  to  500  km)  an  additional  spatial  ( isoplanaticsm)  delay 
exists.  Filter  estimates  are  for  turbulence  effects  along 
the  line-of-sight ; if  a significant  aim-ahead  is  required  to 
fire  at  a moving  target,  the  beam  will  be  aberrated  by  a sec- 
tion of  the  atmosphere  displaced  from  the  line-of-sight  by  the 
lead  angle.  A predictor  similar  to  that  used  to  compensate 
transit  time  delays  in  the  active  multidither  estimator  could 
be  implemented  to  account  for  the  aim-ahead  delay.  Using  a 
r-  spectral  description  of  the  spatial  characteristics  of  the 
turbulence,  a Markov  model  approximation  could  constructed; 
using  this  Markov  model,  the  line-of-sight  turbulence  estimates 
could  be  propagated  (predicted)  through  the  lead  angle. 

Disturbance  Model  Refinement  - The  success  of  the  op- 
timal estimation  and  control  approach  is  attributable  to  the 
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use  of  a precise  mathematical  definition  of  the  disturbance 
dynamics.  In  the  present  study,  the  turbulence  disturbances 
were  characterized  as  five  independent  random  processes.  In 
actual  fact,  the  modal  turbulence  processes  are  cross-correlated; 
this  is  reasonable  since  the  five  turbulence  aberration  modes 
are  all  produced  by  the  same  atmosphere.  Knowledge  of  these 
cross-correlations,  based  on  cross-spectra  for  the  five  modes, 
for  example,  could  be  incorporated  into  the  state  estimation 
algorithms;  this  additional  information  could  be  particularly 
advantageous  to  the  active  multidither  estimator  which,  in  the 
present  study,  must  estimate  the  states  of  five  independent 
processes  from  a single  measurement.  Cross-spectra  derived 
from  a mathematical  model  similar  to  that  used  by  Hogge  and 
Butts  (Ref.  1)  would  provide  the  basis  for  a detailed  charac- 
terization of  modal  turbulence  dynamics  including  cross-correla- 
tions of  the  modal  dynamics. 

Phase  Conjugate  System  - Efforts  to  refine  the  phase 
conjugate  system  should  be  directed  towards  the  shearing  inter- 
ferometer. As  Fig.  5.1-1  showed,  reduction  of  sensor  noise 
level  results  in  a commensurate  improvement  in  Strehl  ratio. 
Optimization  of  the  spatial  filtering  of  the  shearing  inter- 
ferometer, i.e,,  optimal  number  and  location  of  subapertures, 
is  another  potential  area  for  system  improvement. 

Development  of  a Practical  Algorithm  - In  a real  time 
implementation,  problems  arise  which  are  not  given  close  atten- 
tion in  the  preliminary  design  stage.  Specifically,  the  algo- 
rithm computation  time  and  the  transport  lag  it  introduces 
becomes  a primary  concern.  In  addition  to  attempting  to  de- 
velop the  most  efficient  software  possible,  a tradeoff  between 
algorithm  complexity  (i.e.,  disturbance  model  order)  and  com- 
putation lag  may  be  necessary  to  "optimize"  a real  time  imple- 
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mentation.  The  effects  of  computation  lag  can  be  artificially 
introduced  into  a non-real- t ime  computer  model,  allowing  sen- 
sitivity studies  to  be  carried  out.  The  true  computation  lag 
must  be  assessed  through  actual  implementation  of  the  algorithm 
on  the  control  system  computer. 


Test  Bed  Implementation  - The  results  of  Section  3 
shown  that  the  Kalman  filter  approach  to  phase  conjugate  es- 
timation yields  a considerable  performance  improvement  over 
classical  estimation  methods.  Because  of  the  superior  perfor- 
mance offered  by  the  optimal  estimation  and  control  approach, 
serious  consideration  should  be  given  to  its  use  in  test  bed 
implementations  (i.e.,  the  adaptive  optics  portion  of  the 
Shared  Aperture  Medium  Range  Tracker). 
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A.l 


LASER/ ADAPTIVE  OPTICS  MODELS 


A. 1.1  Tilt  Mirror 


The  mechanical  model  of  the  tilt  mirror  for  pitch 
motion  is  shown  in  Fig.  A. 1-1.  Table  A. 1-1  lists  the  model 
parameters  for  the  open  loop  dynamics. 


K-2MM2 

VOICE  con 
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MIRROR  STATE  VECTOR: 


0 - MIRROR  SURFACE  ANGLE 

8 

tt  . MOUNT  compliance  angle 
m 

- NORMAM7.cn  VOICE  COIL  TORQUE 


INPUT: 

0„ur  * TILT  CONTROL  COMMAND 
CMU 

DISTURBANCES 

0.  - BASE  MOTION  DERIVED  FROM 
b 

VIBRATION  MODEL 


Figure  A . 1-1 


Tilt  Mirror  Mechanical  Model  and  Variable  Defi- 
nition for  Mirror  Pitch  Motion 


TABLE  A.I-1 

TILT  MIRROR  MODEL  PARAMETERS 


T-1856 


SYMBOL 

DEFINITION 

VALUE 

“h 

Surface-flexure  hinge  natural 
frequency 

(2H)UBz 

Mount  natural  frequency 

( 2 H )oOQQHz 

I ’ 

Mirror  inertia/mount  inertia 

0.5 

T1 

Voice  coil  time  constant 

0 . 0002sec 

^nsax 

Maximum  mirror  displacement 
-|Ds  V 

2mrad 

9b 

Rras  level  of  base  vibration 

(used  0. ) 

input.  Equivalently:  rms 
level  of  mirror  surface  jit- 
ter due  to  laser  vibration 

- 

•I 


The  following  control  objectives  will  be  the  basis 
of  inner  loop  design  for  the  tilt  mirror. 

• Extend  mirror  surface  response  bandwidth 
to  1000Hz  (-6280  rad/sec) 

• Maintain  a stable  mount  compliance 
mode  response. 

A. 1.2  Focus  Mirror 

Figure  A. 1-2  shows  the  mechanical  moael  of  the  move- 
able  secondary  mirror  of  the  focus  mirror.  Table  A. 1-2  lists 
the  model  parameters  for  the  open  loop  dynamics. 

The  following  control  objectives  specify  the  inner 
loop  design  for  the  focus  mirror. 
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mirror  state  vector 


xs  = mirror  surface  displacement 

x - MOUNT  COMP!.!  ANCE  DISPLACEMENT 

m 


NORMALIZED  VOICE  COlt,  FORCE 


INPUT  : 

U - FOCUS  CONTROL  COMMAND 


DISTURBANCES 

x.  - BASE  MOTION  DERIVED  FROM 
VIBRATION  MOPE!. 


Figure  A. 1-2  Focus  Mirror  Mechanical  Model 


TABLE  A. 1-2 


FOCUS  MIRROR  MODEL  PARAMETERS 

T-1857 


symbol 

DEFINITION 



VALUE 

ws 

Movable  mass-spring  natural 
frequency 

92 . 7 rad/ sec 

m 

Mount  natural  frequency 

(27.  ) 2000Hz 

. 

m 

Mirror  inertia/mount  inertia 

0 . 5 

- 

Voice  coil  time  constant 

0 . 0002sec 

d„,  „ 

Maximum  mirror  displacement 

(used  ») 

max 

< X - X 

- 1 s m 

Primary  mirror  focal  length 

Ira 

a 

Rms  level  of  vibration  input. 

( used  0.  ) 

Xb 

Equivalently:  rms  level  of 

mirror  surface  jitter  aue  to 
laser  vibration 

( used  0. ) 


Extend  secondary  mirror  response 
bandwidth  to  500Rz  - (3140  rad/sec) 


• Maintain  a stable  mount  compliance 
mode  response. 

A. 1.3  Deformable  Mirror 

The  equations  describing  the  dynamics  of  the  defor- 
mable mirror  model  developed  in  Section  2.2.3  are  given  by 
Eq.  A. 1-1  to  A. 1-4.  Table  A. 1-3  lists  the  primary  and  coupling 
dynamics  parameters 

z “ F(z)z  + Bu  (A. 1-1) 

z * vector  of  mode  strength  and 
time  derivative  of  each  mode 

u * vector  of  mode  commands 


TABLE  A. 1-3 

DEFORMABLE  MIRROR  PRIMARY  DYNAMICS  PARAMETERS 

T-1839 


SYMBOL 
■ ■ — — . 1 

DEFINITION 

VALUE 

! 

2 ' 
~V  * I 

Ko 

i 

Zero  displacement  natural 
frequency  squared 

(2.U)2(  3300Kz)2 

m 

Nonlinearity  constant 

( used  0 . ) 

dk  * 

Modal  damping  coefficient 

•2*k 

Ko 

Z 1 

max  > 

Maximum  mirror  displacement 

10um 

¥ 


For  the  kth  mode : 


J2k-1 


'2k 


42k 


uik(z2k-l)  2 z2k-l  " dkz2k 


Primary  Dynamics 


n 


+ 1 0kJzJ  + 

j?2k-l,2k 

Coupling 


(A. 1-2) 


“k<z2k-l> 


Excitation 
(A. 1-3) 


wk(z2k-l)  2 * ^kQ  + (u)k}iz2k-l 


(A. 1-4) 


A. 1.4  Thermal  Deformation 


Parameters  required  by  the  thermal  deformation  model 
of  Section  2.2.5  are  given  by  Table  A. 1-4. 

A. 1.5  Laser  Induced  Vibration 

Figure  A. 1-3  shows  the  spectral  representation  of 
the  laser  induced  vibration,  its  second  order  shaping  filter 
approximation,  and  dynamic  system  realization  as  developed 
in  Section  2.2.4. 

Parameters  to  be  used  in  the  shaping  filter  are: 


iun  = (2tt)(32Hz)  (A. 1-5) 

; * .1  (A. 1-6) 

The  variance  of  the  noise  input  to  the  filter  is  taken 
as  unity;  the  DC  gain  of  the  filter  is  also  unity.  The  distur- 
bance "input  gains"  of  the  active  optics  components  are  adjusted 
to  achieve  a desired  rms  jitter;  e.g.  in  Table  A. 1-1. 
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TABLE  A. 1-4 

THERMAL  DEFORMATION  MODEL  PARAMETERS 

T-1861 


SYMBOL 

DESCRIPTION 

■ 

I VALUE 

Jave 

_ 

Average  incident  beam  intensity. 
Derived  from  power  levels  cited 
in  Ref.  5 and  internal  beam 
radius 

1 

i 

rint 

Internal  beam  radius 

0.10  m 

a 

Surface  absorption  coefficient 

, 0.001 

i 

5i 

Beam  incidence  angle 

45°  all  mirrors 
except  beam  ex- 
pander. 90° 
beam  expander 
primary  and 
secondary 

* 

Radiation  wavelength 

3.8  um 

’1  ’ %2 

Mirror  property  dependent  factors. 
Determined  from  empirical  relation- 
ships in  Ref.  3 

( used  0 . ) 

- 

‘‘thermal 

Thermal  response  time 
(cooled  mirrors) 

. 1 sec 

£ 

thermal 

Thermal  growth  rate 
(uncooled  mirrors) 

.003  um/sec 

A. 1-6  Relay  Mirror  Jitter 


The  equation  of  motion  for  vibration  induced  pitching 
of  the  i relay  mirror,  as  developed  in  Section  2.2.7,  is  given 


by  Eq.  A. 1-7.  The  corresponding  yaw  equation  of  motion  is  of 
the  same  format  with  ^ replacing  and  subscript  v replacing 
subscript  p. 


SPECTRAL  REPRESENTATION 


DYNAMIC  SYSTEM  MODEL 


GAUSSIAN 
WHITE  NOISE 


SECOND  ORDER 
FITTER 


52‘2{wnS*w2 


LASER 

VIBRATION 


Figure  A. 1-3  Laser  Induced  Vibration  Spectrum  and 

Model  Approximation 


+ 2; 


w 

p,  p 


, 2q 
+ oj  0 . 

pi 


kivb(t+Ti)  (A. 1-7) 


Table  A.l-f  lists  the  parameters  that  will  be  used  in 
this  model. 


In  order  to  simulate  the  relative  excitation  levels  of 
the  individual  relay  mirrors  as  described  in  Section  2.2.7,  the 
k^'s  will  be  generated  by  a random  draw  (on  the  initial  cycle  only) 
from  the  uniform  probabilty  distribution  function  shown  in 
Fig.  A. 1-4. 


1max 


Figure  A. 1-4  PDF  for  Relay  Mirror  Vibration  Gains 
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TABLE  A. 1-5 

RELAY  MIRROR  JITTER  MODEL  PARAMETERS 

T- 1860 


SYMBOL 

DEFINITION 

VALUE 

i 

» H 

i''  relay  mirror  natural  fre- 

(2.1  )5Hz  ; 

Pi 

quencies  in  pitch 

i » 1 to 

4 

itfa  relay  mirror  damping  ratio 

0.9; 

in  pitch 

i * 1 to 

- 1 

“± 

^v 

i^  relay  mirror  natural  fre- 

(2H)5Hz; 

yi 

quency  in  yaw 

i ■ 1 to 

4 

’v 

* h 

iz  relay  mirror  damping  ratio 

0.9; 

yi 

in  yaw  * .01;  i * 1,4 

i « 1 tO 

4 

k. 

th  1 

i relay  mirror  vibration  input 

* 

X 

gain.  Randomly  chosen  from  a 
uniformly  distributed  probability 
density  function  on  initial  cycle.  j 

j 

' i 

! 

i relay  mirror  vibration  input  ; 
phase  shift  randomly  chosen  from  : 
a uniformly  distributed  probability  J 
density  function  on  initial  cycle.  | 

21 

^Pi 

k.  is  determined  from  the  maximum  expected  mirror 
max 

jitter  variance  by  working  backwards  through  the  dynamics  of 
the  relay  mirror  and  laser  induced  vibration  models. 


A . 1 . 7 Beam  Expander 


Table' A. 1-6  lists  the  parameters  that  will  be  used 
tc  describe  beam  expander  effects. 


tm*  1 ■? 


TABLE  A. 1-6 

BEAM  EXPANDER  PARAMETERS 


i! 

f 


A. 2 ATMOSPHERIC  PROPAGATION  EFFECTS 

A . 2 . 1 Atmospheric  Turbulence  Models 

Table  A. 2-1  lists  the  shaping  filter  parameters  to 
be  used  in  simulating  atmospheric  turbulence  effects  for  the 
AMD  scenarios;  Table  A. 2-2  lists  the  filter  parameters  for 
the  PC  scenario.  When  driven  by  unity  variance  white  noise, 
the  spectral  content  of  a given  filter  output  approximates 
spectra  generated  by  tne  methods  of  Hogge  and  Butts  (Ref.  1). 


TABLE  A. 2-1 

SHAPING  FILTER  PARAMETERS  FOR  AMD 


GENERAL  FORM  OF  TRANSFER  FUNCTION:  K~£  7- ..  w , 

a ( s-^a ) (s+b  ) ( s^-c  ) 


MODE 


X 

4.4* 10" ~ 

1 

! 14.14 

t 

| 659.5 

! 1 

; 659.5  ! 

1 ' 

X 

7 

3.9xl0“7 

14 . 14 

659.5 

659.5  i 

X 

2.2  R2 
x *1  - T I 

7.2*10‘8 

28.25 

; 

1005. 

| 1005.  j 

X 

2 2 

x -y  1 

1 

-7  ' 

1.3x10  | 

. 

12 , 55 

•j 

to 

0 

773.0  i 

x 

xy  | 

1.2  10"9 

1.410 

j 22.00  1 

880.0  1 

1 

. 0157 

TABLE  A. 2-2 

SHAPING  FILTER  PARAMETERS  FOR  PC 

T-1863 

GENERAL  FORM  OF  TRANSFER  FUNCTION : K 


d (s+a)(s+b)(sfc) 


1 

I 

i 

t 

A. 2. 2 Speckle  Model 

Figure  A. 2-1  is  a block  diagram  of  the  speckle  model 
developed  in  Section  2.3.3. 


w 


R-29294 


MEASURED 

RETURN 


IDEAL  TARGET 


RETURN 


Figure  A. 2-1 


Speckle  Dynamic  Model 


For  the  AMD  scenario : 

a - 211  ( 2300Hz) 


aspeckle  ' ^ 


(A. 2-1) 
(A. 2-2) 


A. 2. 3 Receiver  Model 


